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VULNERABILITY REDUC TION USING MOVEMENT AND SHELTER 


CHAPTER I 
INTRODUCTION 

The purpose of this study is indicated by the following quotation from 
the Scope of Work statement in the Contract: 

"Specific work and services to be performed shall 
include, but are not limited to, die following: 

1. Investigate the relationships among vulnerability, 
warning times, shelter costs and other factors involved in 
strategic movement using basic data developed under Contract 
OCD-OS- 62 -248. 

2. Postulate various combinations and mixe^ of stra¬ 
tegic evacuation and shelter and analyze and compare them. 

3. Devise methods cf evaluating overall, alternative 
plans for reducing vulnerability on the basis of survival rates, 
warning times, costs, time required for activation and other 
factors 

4. Examine the evolutionary development, charac¬ 
teristics, and desired order of development of sui rival 
capabilities related to vulnerability reduction. 



As indicated in Par. 1 of the above Quotation, the present study if 
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a continuation of earlier Dikewood strategic movement studies. ' The 
overall effort may be placed in perspective with the aid of Fig. 1 and some 
intuitive arguments. 



Fig. 1 

survivors Versus Time for 
Various Movement and Shelter Policies 


With very long action times (time between 
lethal effects) complete evacuation of a city cou 


decision to act and arr’val 
Id be carried out, and the 


P J. Flanagan, et al. . Specific Strategic Movement St 
Corporation Final Reoort on Contract OCD-OS-62-248. 
May ’963. (Confidential' 


udies , Dikewood 
DC -FR- 1030; 


S D. Stearns, A Mathematical ?.* ->ael ;or Strategic Movement, 
Ope ra tions Research , Voi. 12, No. 2; March-April 1964 
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immediate survival rate would be essentially 100%, if the evacuated city 
v/ere bombed This survival rate must be reduced by fallou casualties in 
the reception area, but given sufficient time, improvised fallout protection 
can be provided, and a very high attack survivor rate can be achieved 
(Curve A in Fig 1). If fallout protection were already available in the recep¬ 
tion area, the action time required to achieve this survival rate would be 
reduced (Curve B) If shelters (from initial effects as well as fallout) were 
provided in the target area and the population vulnerability decreased for 
short action times, there would still be enough people killed either directly 
or indirectly that the total survival rate would probably be smaller than that 
possible with evacuation and a long action time (Curve C) Combinations 
that evacuate those in the most vulnerable locations while providing shelter 
for these m somewhat safer positions might lead to survival curv:s like the 
dashed curves in the figure (Curve D) 

The development of realistic curves of the types shown in Fig 1 
may be done in three major steos 

1 Postulation of alternative movement and shelter policies. 

2 Development f movement and shelter plans based on 
the policies in (1), and 

3 Evaluation of plans developed in (2) against the 
range of attack conditions considered reasonable 

In carrying out the first step one would hope to postulate as many 
policies as imagination allows, this is largely an intuitive step 



Once a policy is postulated, a technique for developing plans to car¬ 
ry out the policy can be constructed. In this report. Chapter JI describes a 
basis for development of strategic movement plans previously described 
in more detail in Refs. 3 and 4. Chapter III describes tools developed to 
fulfill Step 3 for strategic movement plans. 

Chapter IV describes some possible shelter policies (Step 1) and 
computer programs that can be used as a basis for development of shelter 
plans (Step 2). The evaluation of shelter plans (Step 3) against various 
attacks is a fairly straightforward step, but the procedure for doing so 
has not been automated. 

Chapters II-IV then are aimed at the final goal of developing pre¬ 
ferred mixtures of movement and shelter typified by an evaluation curve 
of Type D in Fig. 1. Another independent approach to the development 
of such preferred mixtures was also followed and is described in Chapter V 
and in Ref. 5. 

All four items in the Scope of Work statement of this contract 
describe problems that will require continuing treatment by the Office of 
Civil Defense. The effort summarized in this report provides some help 

3 

R J. Flanagan, et al., Large-Scale Strategic Movement Planning , 
Dikewood Corporation Technical Note DC-TN-1039-1; January 15, 1964. 

4 

S. H. Dike, et al. , A Computer Program for Planning a Strategic 
Movement , Dikewood Corporation Technical Note DC-TN-1039-9; 

May 24, 1965. 

5 

K. D. Granzow, A Model for Development of Preferred Mixtures of Evac¬ 
uation and Shelter, Dikewood Corporation Technical Note DC-TN-1039-2; 
July 6, 1964, 
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in understanding Pars. 1, 2, and 3, but very little toward Par. 4. Some 
ideas for treatment of problems associated with Par. 4 are being developed 
and will be reported later under Contract OCD-PS-65-53, where further 
work in these areas is being supported. 
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CHAPTER II 


LARGE-SCALE STRATEGIC MOVEMENT PLANNING 

1.0 Introduction 

Capability for movement has probably always been understood to 

be a proper element of a defense system. Its potential value in many 

kinds of thermonuclear conflicts has been pointed out by a few workers, 

notably Herman Kahn and some of his former co-workers at RAND. 

The Hudson Institute staff examined some of the variables involved 

6 

in strategic movement problems as a part of an OCD-sponsored study. 

In that study, a set of three illustrative plans was prepared for the north¬ 
eastern section of the U. S These plans were associated with various 
levels of crises, the primary effort being devoted to an evacuation that 
would take a week to complete. Modifications to the basic one-week plan 
to illustrate some effects of other warning times were also examined 
These included a two-day plan and a plan that might be associated with a 
crisis that escalates over a one-month period This set of plans provided 
insight into a number of problems anu has proved very valuable in Dike- 
wood studies 

The Dikewood Corporation was asked to perform a study of strategic 
movement from two cities: one city that is relatively isolated and one city 

^William M Brown, Editor, Strategic and Tactical Aspects of Civil 
Defense with Special Emphasis on Crisis Situations , Hudson Institute 
Report No HI-160-RR, January 7, 1963 
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surrounded by other population centers so that evacuees from it would have 
to compete with other evacuees for reception area space The two cities 
chosen were Albuquerque and New York The development of plans for 
these two particular cities included an examination of some questions con¬ 
cerning the feasibility of strategic movement. Some of the factors con¬ 
sidered in Ref 1 were re-examinea and are discussed in this chapter. 

In the course of the study, a computer program was developed that 
can be used t.o calculate the evacuee distribution for which expected casu¬ 
alties are a minimum within an attack area that includes a number of 
evacuation and reception sectors. This program is referred to as the 
Distribution program and is described briefly in this chapter and in 
more detail, in Refs 1, 2, and 3 

Two alternative techmqu- s for planning movements to achieve any 
desired distribution have been developed and are also described in this 
chapter The planning techniques make use of either of two computer 
programs referred to as the Movement and Min Man-Mile programs The 
Distribution program can be used to decide how many evacuees should be 
housed in each reception place and the Movement and Min Man-Mile pro¬ 
grams can be used to decide where they come from when they should 
leave, and what routes are to be followed 

The techniques discussed in this chapter are not suggested as being 
unique, optimum or in any sense the only right way to plan strategic move¬ 
ment However, they do represent one method that has received much 



careful consideration and is as nearly automated as seems reasonable. 
Again, computers are not expected to add any magic to the results, but they 
do provide an accurate, convenient, and inexpensive means of handling the 
large amount, of bookkeeping required in planning such an operation. 

If the planning techniques described in this chapter were applied, the 
bare minimum of planning needed tor an emergency capability to perform 
a strategic movement would be established. However, much of the detailed 
planning required to make one confident that the operation would be suc¬ 
cessful and that it would proceed smoothly would require further consider¬ 
ation. These details should be treated later by planning groups working, 
for example, at the state or local level. These more detailed plans might 
be prepared in the same way that the State Survival Plans were prepared, 
in the way the National Fallout Shelter Survey was performed, or in some 
similar manner. 

2.0 Prediction of Expected Casualties and Their Relation to the 
Evacuee Housing Problem 

One reasonable criterion for use in planning a strategic movement 

is the minimization of the expected number of casualties, and a previous 

Dikewood effort emphasized the development of a technique for calculating 

1 2 

such minimum-casualty distributions. ' To find a minimum-casualty 
distribution for any given area such as a state or group of states, the area 
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was divided into constant-vulnerability sectors in which the initial population 
and shelter protection factors and capacities were known. The number of 
people that can move from one sector to all others in the time period of 
interest was then specified. With these data items it is possible to calculate 
the minimum-casualty distribution. A computer program was developed to 

facilitate the calculation; the program is described briefly in the final report 

1 2 
on Contract OCD-OS-62-248 and in more detail in an earlier report writ¬ 
ten under the same contract. The technique was applied to Albuquerque and 
New York City to develop specific distribution plans and to see how the 
choice of plans affects the numbers of casualties. Albuquerque was treated 
as an isolated city, i. e., a particular reception area was chosen arbitrarily 
and various plans for moving people within this one area were studied. For 
New York City, the problems associated with competition for space were 
emphasized by examining plans for various groups of states in the Northeast. 

In the studies of both of these areas, fatalities were related to housing 
load factor. Housing load factor is a measure of billeting burden and is 
defined as the ratio of population after movement to the population prior to 
movement in any particulai area of interest. For example, if some town 
experiences a load factor (LF) of two, it simply means that people are 
"doubling up. 

Calculations of fatalities verified the intuitive expectation of a large 
decrease in fatalities associated with emptying pl-ice c *hat receive ini^ai 




effects. For example, the relation between load factor and fatalities asso- 

7 

ciated with the NAHICUS '63 attack is shown in Table I for Albuquerque 
and its reception area . 


Table I 

Fatalities versus Maximum Housing Load Factor for Residents 
}f Albuquerque and Its Reception Area 


Maximum housing 
load factor 

Expected 
fatalities (%) 

First 

differences 

1 (no evacuation) 

61 


2. 73 (uniform) 

12 

28 

3 

10 

7.4 

5 

8 

1 

7 

7 

0. 5 

9 

7 

0 


a 


For a unit change in maximum housing load factor. 


There are several reasons for preferring a load factor near uniform. 
First, the results in Table I indicate that there is still a sizeable gain for 
increasing tht load factor to a level somewhat greater than uniform, but 


7 

Nuclear Attack Hazard in Continental U. S. , 1963 , Office of Emergency 
Planning and Department of Defense. 

a. Vol. I, Problem and Approach (Confidential). 

b. Voi. II, Methodology and Input (TSRD). 

c. Vol. Ill, Summary of Attack-Effects Probabilities (Secret). 
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that tlie rate of reuirn then decreases rapidly (See first differences in 
Table I ) This is not unexpected because increasing the load factor above 
the uniform level simply means that the areas expected to experience 
heaviest fallout will be avoided Also, strain upon both hosts and evacuees 
decreases with decreasing load factor Finally, three important sources 
of uncertainty in the calculations upon which Table I is based also lead one 
to favor load factors near uniform. These uncertainties are 
a Fallout prediction techniques, 
b Winds, and 
c. Enemy intentions. 

It is fairly well-known that various techniques for calculation of 
fallout patterns even for a single wind structure may lead to some rather 
large differences in predicted dose patterns. This is at least partly be¬ 
cause there has never been conducted a well-instrumented large-yield 
burst over a land surface with strong winds aloft For obvious reasons, 
such a test may never be carried out even though it would do much to clarify 
the fallout ed< non problem 

Unce rtaintles m wind structure are another source of difficulty 

One example that illustrates this point is given m Fig 2 where the 1000-r 

dose contours for winds observed on four particular days are shown or the 
8 

same map Figure 2 indicates that, while it may be reasonable to plan on 

g 

DAY oun g- lout , <>. - at Albuquerque. New Mexico Sandia 

Corpora?ion Technical Memorar. ’j m SCTM-1 95 - 59(51). January, 1 960 





Note \Veapon yield = 2C Mt, .alf fission 



Fig 2 

San ia Fallout Model Predictions of the 100C-R Dose 
Contours for Bursts at Albuquerque on Four Particular Davs 
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expected value basiu, one should be aware that almost any wind pattern is 
at least possible. An investigation of fallout wind statistics may be found 
in Refs. 9 and 10. 

The third source of uncertainty mentioned above, the enemy's 
intentions, has been discussed at length m a number of publications (e g , 
see Ref. 11). As a hedge against the particular uncertainty of whether an 
enemy would consider population a principal target, a bonus target or a 
target to be spared, one might employ a plan that distributes the popu¬ 
lation in a nearly uniform manner. For a fixed hardness posture such a 
distribution would tend to impose the greatest cost upon an enemy that has 
as a goal the destruction of some large fs cuon of the T\ S. population. A 
uniform distribution can be approached if some upper limit is placed on the 
size of post-movement population groups 

Next, it is appropriate to consider the numerical values ot uniform 
load factor for a large scale evacuation and to consider whether they ar^ 
prohibitively iargt One approximation can be made by examining the frac¬ 
tion of the l. S population contained ;n cities versus i tty size Then, if 


E D Callahan, et al. . The Frobable Fallout Threa’ Over the Conti¬ 
nental United States Technical Operations, Inc. Report No TO-B-60-13. 
Dece mber i, 1960. 

8. N Charles, M ean Lay er Wi nds by Seasons . Sand;a Corporation and 
t. . S. Weather Bureau Cooperative Project in Climatology Phase Two, 
March 1960. (unpublished) 


Hermar Kahn. Tmnk-r.g about th° Cr.thmkable , Horizon Press N V . 1962 


Strategic evaluation allows 
is a.med at the destruction 
would also k:’l a large frac 


the offense plan a countervalue vtac* ‘hat 
of fac.l;t;es without requiring an attack that 
tior. of the population 
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one can decide on the minimum size of a city to be evacuated, an estimate 
of the uniform load factor can be made on a national basis. For this esti¬ 
mate the Census Bureau's "urbanized areas" have been used because they 
are believed to provide the best separation of urban and rural population 
in the vicinity of large cities. The results are shown in Table II. Note 
that even if all urbanized areas with populations greater than 50, 000 were 
evacuated, the uniform load factor is not much more than 2. 


Table II 

Estimate of National Uniform Load Factor 
for Evacuation of U. S. Urbanized Areas'* 

Minimum population Fraction of U. S. Uniform 

— (thousands) population load factor 


1000 

. 29 

1. 4 

7 50 

. 33 

1. 5 

500 

. 37 

1. 6 

2 50 

. 43 

1. 8 

2 00 

. 46 

1. 9 

ICO 

. 51 

c. 0 

75 

. d2 

2. 1 

50 

. 54 

0. 2 


a Data from I960 U S Census. 

Of course, this is a very rough approximation, but it serves to indicate that 
me problem is soluble A more detailed examination of the problem has 
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been made and it tends to bear out the same conclusion. In this second 
look, the load factor was calculated for each state, since it is believed 
that the smaller the number of gove r nments involved, the more manage¬ 
able the plan. Thus, it. seems preferable ;o have each state stand alone or 
where this is impossible, to have the smallest reasonable group of states 
involved in supporting each other. While this policy of minimizing the num 
bet of interstate agreements required may cause some hardship, it is 
expected that it would lead to a workable plan in the shortest time. 

As a further improvement on the first load-factor estimates, a 
specific heavy attack was chosen and the load factors were calculated 
assuming that everyone !’v:ng within 24 miles of each target would be 
evacuated. This radius was originally chosen somewhat arbitrarily on 
the basis that it is the distance to an overpressure of 1 psi plus 3 times 
an estimated CEP 1 5 miles for a surface burst of a 20-Mt weapon. 
These choices were expected to lead to a conservative estimate of the 
number of evacuees since ♦hev associated with a 0. 002 probability of 
exceeding 1 psi on the circumferenct. The results of the analysis of uni¬ 
form ioad factors by spates for the Te< h Ops attack are shown m Fig 3. 
Note that some of the states will have to be grouped, panuu’urlv in the 
Northeas*. 

The chosen attack is often referred to as the Tech Ops Attack" and 
includes dehvfry of 816 weapons on l\ S. rrr’.jry, industrial, and 
da.r. tnrgets (Hef d) 
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Fig. 3 

Uniform Load Factor fo^ Each State 


-91- 















Because the results of the analysis depicted in Fig. 3 seemed favor¬ 
able, a post-movement limit of 1 00, 000 people was imposed on reception 
city size and the uniform load factors for the remaining areas were recal¬ 
culated. This calculation was performed for each state and for several 
groups of states in the Northeast. The results are given in Tables III and 
IV. It should be pointed out that no check was made to be sure that there 
were no more than 100, 000 people in groups located within 24 miles of each 
other. 

Some indication of the size of load factors that might become accept¬ 
able under the desperate conditions associated with a "total" evacuation can 
be obtained by comparing present U. S. housing levels to those of other 
countries. Such a comparison is made in Table V where it may be noted 

that current load factors of 3 are common and that they are as high as 5. 2 

12 

relative to U. S. standards. One might guess that load factors larger 
than twice the maximum current world value would be tolerable for only 
a short time or under very desperate conditions. The data in Tables III 
and IV indicate that a maximum load factor of less than ten is eaoiiy 
achieved by grouping states and that even a maximum load factor of less 
than five may be achievable everywhere in the U. S. 

12 ” 

United Nations Statistical Yearbook for 1961 , Statistical Office of the 
United Nations, Department of Economic and Social Affairs, 1961. 


Table III 


Population and Load Factor Data for Conterminous U. 


S. 


a 



Total b 
population 

Target 

area 

population 

Reception 

area 

c 

population 

Q 

Uniform load factor 
No limit*"* 10 5 limit 

Alabama 

3,267 

1,218 

2, 049 

1.6 

1.6 

Arizona 

1, 302 

794 

508 

2. 6 

2. 6 

Arkansas 

1, 786 

451 

1, 335 

1.3 

1.3 

California 

15, 717 

12,088 

3, 629 

4. 3 

4. 6 

Colorado 

1, 754 

551 

1,203 

1. 5 

1.5 

Connecticut 

2, 535 

2,470 

65 

39 

39 

Delaware 

446 

377 

69 

6.4 

6.4 

District of 
Columbia 

764 

764 

0 

- 

- 

Florida 

4, 952 

3, 180 

1, 772 

2.8 

2.9 

Georgia 

3, 943 

1, 644 

2, 299 

1. 7 

1.7 

Idaho 

233 

12 

221 

1. 1 

1.1 

Illinois 

10, 081 

7, 136 

2, 945 

3. 5 

3. 6 

Indiana 

4, 662 

1,978 

2, 684 

1. 7 

1.8 

Iowa 

2, 758 

578 

2,171 

1. 3 

1.3 

Kansas 

2, 179 

1, 136 

1,043 

2. 1 

2. 1 

Kentucky 

3, 038 

1,074 

1, 964 

1.5 

1.5 

Louisiana 

3, 257 

1, 767 

1, 490 

2.2 

2.2 

Maine 

969 

680 

289 

3.4 

3. 4 

Maryland 

3, 100 

1, 504 

1, 596 

1.9 

1.9 

Massachusetts 5, 148 

5, 115 

33 

15* 

- 

Michigan 

7, 823 

5, 853 

1,970 

4.0 

4.2 

Minnesota 

3, 414 

1, 161 

2,253 

1.5 

1.5 

Mississippi 

2, 178 

123 

2,055 

1. 1 

1.1 
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Table III (Continued) 


cl 

Population and Load Factor Data for Conterminous U. S. 



Total b 
population 

Target 
area ^ 
population 

Reception 

area 

c 

population 

c 

Uniform load factor 

No limit d 10 5 limit d 

Missouri 

4, 320 

2,467 

1, 853 

2. 3 

2.5 

Montana 

675 

112 

563 

1.2 

1.2 

Nebraska 

1,411 

597 

814 

1. 7 

1.7 

Nevada 

285 

79 

206 

1.4 

1.4 

New 

Hampshire 

607 

272 

335 

1.8 

1.6 

New Jersey 

6, 067 

5, 684 

383 

16 

19. 3 

New Mexico 

951 

468 

483 

2. 0 

2.0 

New York 

16, 782 

14,977 

1, 805 

9. 3 

10. 6 

North 

Carolina 

4, 556 

1,252 

3, 304 

1.4 

1.4 

North Dakota 

214 

103 

111 

1.9 

1.9 

Ohio 

9, 706 

7, 393 

2, 313 

4. 2 

4. 5 

Oklahoma 

2, 328 

966 

i, 362 

1. 7 

1.7 

Oregon 

1, 769 

754 

1,015 

1. 7 

1. 7 

Pennsylvania 

11, 319 

8, 440 

2, 879 

3. 9 

4. 1 

Rhode Island 

859 

859 

0 

- 

- 

South 

Carolina 

2, 383 

930 

1, 453 

1. 6 

1. 7 

South Dakota 

681 

107 

574 

1.2 

1.2 

Tennessee 

3, 567 

1, 998 

1, 569 

2. 3 

2.3 

Texas 

9, 580 

4, 868 

4, 712 

2.0 

2. 1 

Utah 

891 

576 

315 

2 8 

2.8 

Vermont 

390 

95 

295 

1. 3 

1.3 

Virginia 

3, 967 

1. 822 

2, 145 

1.8 

1.9 
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Table III (Continued) 


Population and Load Factor Jata for Conterminous U. 


S. 


a 



Total k 
population 

Target 

area 

population 

Washington 

2,853 

2. 016 

West Virginia 1,860 

510 

Wisconsin 

3,952 

1, 671 

Wyoming 

330 

54 

Totals 

177,609 

110, 733 


Reception 

c 

Uniform load factor 

area 

c 

population 

No limit** 

10^ limit** 

837 

3.4 

3.5 

1, 350 

1.4 

1.4 

2,281 

1. 7 

1. 7 

276 

1.2 

1.2 


66, 876 2.7 


a 

Assumes evacuation of all people within a 24-mile radius of the targets in 
the Tech Ops list of military, indus + rial, and dam targets (Ref. 9). 

k All populations are expressed in thousands. 

Uniform load factor equals "Total population" divided by "Reception area 
population. " 

^Refers to post-movement limit on size of places in reception areas. 
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Table IV 


Housing Load Factors for Northeastern States 


Uniform load factor 



Group 

No Lmit a 

10 5 limit' 

I. 

New York, New Jersey, 

Pennsylvania 

6. 7 

7.2 

II. 

Group I plus West Virginia 
and Ohio 

5.2 

5. 5 

III. 

New York, New Jersey, 

Pennsylvania, West Virginia, 

Virginia, Maryland, Delaware, 
and District of Columbia 

4. 3 

4. 5 

IV. 

Maine, Vermont, New Hampshire, 
Connecticut, Massachusetts, and 
Rhode Island 

10. 3 

11.1 

V. 

Group III plus group IV 

4. 9 

5. 1 

VI. 

Group V plus Indiana, Michigan, 

Ohio, Kentucky, Tennessee, and 
North Carolina 

3. 5 

3. 6 


a 


Refers to post-movement limit on population of places in reception areas. 
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Table V 


Average Housing Load for Various Places 


a 


Plaee a 

Data for year 

No. of 

persons per room 

Load 

factor 

Argentina 

1947 

2.2 

3. 7 

Bulgaria 

1956 

1.8 

3. 0 

Canada 

1951 

C. l 

1.2 

Czechoslovakia 

1950 

1.5 

2. 5 

Denmark 

1955 

0. 7 

1.2 

Dominican Republic 

1955 

1.7 

2.8 

Finland 

1950 

1.5 

2. 5 

France 

1954 

1.0 

1. 7 

Germany, Federal 
Republic 

1956 

1.0 

1. 7 

Greece 

1951 

1.8 

3.0 

Guatemala 

1949 

3. 1 

5. 2 

Italy 

1951 

1.3 

2. 2 

Poland 

1950 

1.4 

2. 3 

Spain 

1950 

1 . 1 

1.8 

USSR 

1956 

1.5 

2. 5 

UK 

1951 

0.8 

1. 3 

Yugoslavia 

1954 

2. 3 

3. 8 

US 

1960 

0. 6 

1.0 


a Source is United Nations Statistical Yearbook, 1961 (Ref. 1?) 
k Load factor is measured relative to the U. S. for 1960. 
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To summarize the remarks on choice of load factor: 


a. The facts that expected casualties decrease slowly with 
load factor above the uniform level, tha* the strain among hosts and evac¬ 
uees decreases with decreasing load factor, and that load factors close to 
the uniform level are a hedge against ui certainties, make a load factor 
slightly greater than uniform the preferred choice. 

b. Load factors smaller than those commonly accepted for 
everyday living by other people in the world are achievable within most 
states or at worst by grouping several s:ates. 

3. 0 Effects of Providing Simple Fallout Shelter 

The numbers of fatalities shown ;n Table I for ♦he Albuquerque area 

were obtained under the assumption that anyone who lacked space in a 

National Fallout Shelter Survey (NFSS) Fhase 1 shelter or a Phase 2 mine 

$ 

shelter was protected only by a house (assumed PF=2). Since there are not 

many shelters in the reception area, this meant that most people were 

protected with a PF of only two. There a .re at least two simple types of 

temporary shelters that can be constructed rapidly to obtain a PF of at 

least 20 and, for a small adcaional cost, a PF of 100. The Tech Ops 

13 

improve^ basement shelter is one type that can be used in many parts 

At 

There are practically no basements in New Mexico. 

13 

E D. Callahan, L Rosenblum, J R Coombe, Shelter from Fallout , 
Tech Operations, Inc., Report No. TO-B60-30, April 1961. 
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of the U. S.; the other is the simple trench. A few rudimentary experiments 
with the construction of trench shelters were performed by two of the authors 
and details which indicate their feasibility are reported in Ref. l;the Hudson 


Institute and Research Triangle Institute have also considered use of trench 


shelters and believe it to be feasible 


6, 14 


Table VI shows the results of calculations of the numbers of fatalities 
if a PF=20 or a PF=100 shelter is provided for all Albuquerque evacuees and 
their hosts. Note that if shelters having a PF=100 were provided, no fatalities 
would be expected. Typical results for some northeastern states are given in 
Table VII. 


4. 0 Transportation Problems 

A feeling for the mobility of the U. S. population can be obtained by 
examining the data in Table VIII, where it is shown that the states all have 
an average of four or less people per automobileZ^Of course, the average 
will be large; in cities that have extensive public transportation systems. 

The five boroughs of New York City have the largest average with 5.9 persons 
per automobile. Boston is next with 4. 5, Philadelphia has 4. 1, and all 

g 

others are 4. 0 or less. In addition, there were about 12 million trucks, 
buses, and publicly-owned vehicles in 1960, probably enough capacity to 


K. E Willis, E. R. Brooks, L. J. Dow, Final Report Crash Civil 
Defense Program Study , Research Triangle Institute, Operations 
Research Division, April 30, 1963. 


Highway Statistics, 1951 , l*. S. Dept, of Commerc e, Bureau of Public 
Roads, U. S. Government Printing Office, 1963. ($1.00) 
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Table V I 


Fatalities for Albuquerque Reception Area 
When Simple Fallout Shelters a-e Provided 

Maximum housing_ Fatalities(%) _ 

load factor Existing shelter Min PF-20 Mm PF -100 


2 7 (uniform) 

12 

CL 

— 

3 

10 

1 5 

0 

5 

6 

0 7 

0 

i 

T 

1 

0 5 

0 

9 

** 

1 

0 5 

0 


a 


Not calculated 


Table VII 


Fatalities in Northeastern States' 


Load Factor 


5 2 (uniform) 


8 

10 


Fa’ali* ms f ' r :> 


Mm PF-20 


10 

8 

8 


Mm PF 100 


* D 


Group for these mlculmiors cor -ists of New York New Jersey Permsy 
van a. Wes’ V irginia, and Ohio 


Not calculated 






Table VIII 


Number of Registered Automobiles by State 


State 

Total 

population 

(thousands) 

Registered 

autos 

(thousands) 

Registered 

trucks 

(thousands) 

People 

per 

a 

auto 

People 

per 

truck 

Alabama 

3, 267 

1,065 

239 

3. 1 

13. 7 

Arizona 

1, 302 

520 

142 

2 5 

9. 2 

Arkansas 

1, 786 

516 

211 

3. 5 

8. 5 

California 

15, 717 

6, 892 

1, 187 

2. 3 

13. 2 

Colorado 

1, 754 

751 

217 

2. 3 

8. 1 

Connecticut 

2, 535 

1,011 

128 

2. 5 

19. 8 

Delaware 

44« 

147 

51 

3. 0 

8. 7 

District of Columbia 

7 64 

486 

20 

4. 1 

38. 2 

Florida 

4, 952 

2, 125 

317 

2. 3 

1 5. 6 

Georgia 

3. 943 

1, 261 

291 

3. 1 

1 3. 5 

Idaho 

233 

264 

121 

0. 9 

1.9 

Illinois 

10,081 

3. 389 

456 

3.0 

22. 1 

Indiana 

4, 662 

1. 719 

363 

2 7 

12.8 

Iowa 

2, (58 

1.089 

2 58 

2 : 

10. 7 

Kansas 

2, 179 

897 

291 

2 4 

7. 5 

Kentucky 

3, 038 

977 

253 

3. 1 

12. 0 

Louisiana 

3,257 

°52 

231 

3. 4 

14. 1 

Maine 

969 

310 

73 

3. 1 

13, 3 

Maryland 

3, 100 

1.034 

144 

3. 0 

21. 5 

Massachusetts 

5, 148 

1, 660 

1 94 

3. . 

26 3 

Michigan 

7, 82 3 

2, 938 

406 

2 7 

19.3 

Minnesota 

3,414 

1. 321 

282 

2. c 

12.1 

Mississippi 

2, 178 

548 

1 92 

4.0 

11.3 

Missouri 

4, 320 

t OQ ? 

* , *v * 

344 

3. 3 

12. 6 

Montana 

675 

266 

i 0 

_ W i 

2. 5 

5. 6 

Nebraska 

1.411 

571 

184 

2 5 

i , 
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Table VIII (Continued) 


Number of Registered Automobiles by State 


State 

Total 

population 

(thousands) 

Registered 

autos 

(thousands) 

Registered 

trucks 

(thousands) 

People 

per 
. a 
auto 

People 

per 

truck 

Nevada 

285 

i41 

45 

2.0 

6. 3 

New Hampshire 

607 

229 

46 

2. 7 

13. 2 

New Jersey 

6, 067 

2, 248 

278 

2.7 

21.8 

New Mexico 

951 

327 

115 

2. 9 

8. 3 

New York 

16, 782 

4, 630 

549 

3. 6 

30. 6 

North Carolina 

4, 556 

1, 455 

339 

3. 2 

13. 4 

North Dakota 

214 

234 

114 

0. 9 

1.9 

Ohio 

9, 706 

707 

444 

2. 6 

21. 9 

Oklahoma 

2, 328 

904 

312 

2. 6 

7. 5 

Oregon 

1, 7 69 

763 

176 

2. 3 

10. 1 

Pennsylvania 

11,319 

3, 805 

552 

3. 0 

20. 5 

Rhode Island 

859 

309 

38 

2.8 

22. 6 

South Carolina 

2 . 383 

737 

155 

3.2 

15. 4 

South Dakota 

681 

2 62 

102 

z 6 

6. 7 

Tennessee 

5, 567 

1, 119 

238 

3. 2 

15. 0 

Texas 

9, 580 

3, 611 

938 

2 7 

10. 2 

Utan 

891 

34 5 

90 

2.6 

9.9 

Vermont 

390 

124 

30 

3 1 

13. 0 

Virginia 

3, 9 67 

1,247 

228 

3. 2 

17. 4 

Washington 

2,853 

1,135 

265 

2. 5 

10. 8 

West Virginia 

1, 860 

490 

125 

3. 8 

14. 9 

Wisconsin 

3. 9 Q 2 

1, 355 

2 75 

2.9 

1 4 4 

Wyoming 

330 

142 

67 

2 3 

4 9 

The numbers of automobiles are for the year 1961, toe popular 

ior. for . 

9 u 0 

Because ot this time difference between the two sources of da T i 

i the hr 

. t*- i 


values are slightiy optimistic. The source of motor vehicle data . - R. • 
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carry out an evacuation solely by this means since it would imply an average 
load of less than 10 people per truck. 

In addition to the U. S. capacity for movement of people by motor 
vehicle, there is a large railroad capacity. In the Hudson Institute plans 
for evacuation of the northeastern states, about 20 percent of the evacuees 
moved by rail. This movement would be primarily by freight car, with 
each car holding 65 people and each train including 100 cars. Movement 
by rail has a number of advantages, a principal one being that it is Les 1 
weather-dependent than movement by automobile. 

5. 0 Recent Examples of Large-Scale Movements 

Both experiment and history indicate that strategic movement is 
feasible. Some of the major attempts to test evacuation techniques are 
described in Refs. 16 and 17. In Operation Rideout, a test evacuation of 
Bremerton, Washington, about 2000 vehicles evacuated the downtown area 
and passed the cwty limits in a half hour. The average speed of the traffic 
columns was 30 mph. 

Operation Green Light was a test evacuation of a 1000-block area of 
downtown Portland, Oregon. In 34 minutes about 29, 000 vehicles and 101,000 
people had left the area; this included 11,000 people who walked. 

16 

Operations Walkout, Rideout, and Sc at, National Academy of Sciences, 
National Research Council, 1955 (unpublished). 

17 


Operation Green Light , Disaster Relief and Civil Defense Office, 
Portland, Oregon, September 1955. 


Hurricane Carla provided an example of the feasibility and possible 

value of employing strategic evacuation. In the evacuation of the Texas 

and Louisiana coastal areas, about 500, 000 people moved distances of as 

much as 350 miles in less than two days. Decisions to evacuate were no* 

made simultaneously over the rather large area evacuated and a check of 

individual city or county movement times indicates that the movement rates 

for Operations Rideout and Green Light were not atypical. For example, it 

took about 6 hours to evacuate 108, 000 people from Jefferson County, Texas. 

The Carla evacuations were performed using normal traffic procedures or, 

in some places, by using all but one lane which was left open for emergency 

traffic. Perhaps the best endorsement for strategic evacuation comes from 

the people and officials involved. The following two paragraphs are typical 

18 

expressions of the people directly involved in the evacuation: 

"The success of the Carla operation left coastal offi¬ 
cials without exception sold on evacuation as a practical, 
cheap, and life-saving device. Agreement was unanimous 
among state and local officials that, if they had listened 
to defeatists and critics of evacuation, thousands of lives 
would have been lost. The Port Arthur CD director said, 

'Anyone who says now that total evacuation is impossible 
is crazy. It was proven, we did it. ' 

"The extent of the success startled even those traffic 
experts who had engineered the operation. The State Direc¬ 
tor of the Texas Department of Public Safety said. If 
someone had tuld me that we could have evacuated between 


18 

Mattie E. Treadwell, Hurricam Carla , Department of Defense, Office of 
Civil Defense, Region 5, U. S. Government Printing Office, December 19C1. 
($0. 55). 
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half a million and 750, 000 people, under the stress we 
had, and not have one fatality or injury, I wouldn't have 
believed him. If someone had told me there'd be no 
panic, I wouldn't have believed him. '" 


6. 0 Elements of a Strategic Movement Plan 

The following are believed to be the basic elements that require con¬ 
sideration in a strategic movement plan: 

a. Decision as to who moves and when. 

b. Delineation of evacuation and reception places. 

c. Transportation, including method, traffic control, refueling, 
treatment of breakdown, and similar details. 

d. Billeting, feeding, and medical care. 

e. Fallout shelter and radiological monitoring in reception areas. 

f Supply. 

g. Command and control. 

h. Communications. 

The choice of what groups of people should move and when should be 
made at a high level, and the announcement should probably come from the 
President. Some local officials may anticipate such an announcement, but 
a national decision would still seem appropriate. One interesting scenario 
in which a sequence of events leads to a strategic evacuation may be found 
in Ref 6 (p. Y-B-l ff. ). Another fictional, but plausible, sequence in which 
there was time to employ strategic evacuation but no capability for it, is 
provided m Pat F-ank's novel Alas, Babylon . 
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The choice of evacuation and reception places and their association 
is the main subject of the planning study suggested in Ref. 3. It is believed 
that making this set of choices is a basic step in the development of a capa¬ 
bility for strategic evacuation. Improvisation might take care of most 
problems associated with strategic evacuation, but the relative simplicity 
of assigning reception areas for each evacuation area makes it seem ex¬ 
tremely unwise to risk the possible mass confusion associated with 
improvisation in the assignment of reception areas to evacuees. 

A gross treatment of transportation problems is described in Ref. 3. 
That is, it is suggested that major transportation routes between evacuation 
and reception areas be catalogued and assigned in an efficient manner, but 
that details of keeping these routes full and the traffic controlled, be plan¬ 
ned with local participants. 

Billeting, feeding, and medical care are believed to be matters that 

should be handled locally. Actually, these functions were delegated by the 

President to the Department of Health, Education, and Welfare except for 

some aspects of food management that were assigned to the Department of 
19 

Agriculture. ^ However, it would seem that OCD at least has a responsibility 
to see that these departments are apprised of strategic movement plans and 
that they then make suitable arrangements to fulfill their assigned roles. 

19 

Executive Order 10958, As Amended; August 14, 1961 and Executive 
Order 11001; February 16, 1962. 


-31- 


Fallout shelter in reception areas can be provided for everyone at 
a very modest cost by using improved Kacerre: ' . .cnch shelters. This 

would require a certain amount of preplanning to determine how many 
shelter spaces are needed in each reception area, how many could be con¬ 
structed with available equipment and materials, and how much additional 
stockpiling is required. The study outlined in Ref. 3 would make a start 
in this direction by indicating how many NFSS spaces would be used, how 
many basements, and how many trenches are required lor the remainder. 

Supply problems and command and control problems can be consid¬ 
ered in the following phases: 

a. Pre-movement, 

b. Movement, 

c Post-movement, pre-attack, 

d. Post-attack, in-shelter, or 

e. Return. 

Both supply and command and control problems will require local- 
level planning. However, specification of the distribution of people requiring 
supplies during the post-movement, pre-attack and the post-attack, in-shelter 
phases would be a necessary input This distribution would provide a basis 
for planning the rerouting of normal supply lines and the stockpile locations 
for supplies for these two phases. 

Communications and radiological momtoring needs are «* cpected to 
be fulfilled by meeting the requirements associated with other parts of the 
national CD capability. Of course this is only a judgment and eventually a 
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study would be required to define these particular needs and to determine 
the feasibility of adapting available systems to this application At the very 
least, some paper planning will be required to make use of the existing and 
planned facilities. 

The 1959 State Survival Plans represent an approximation to the type 
of local planning that is needed. It is believed that sucn local plans can be 
made as they were before, but that procedures should be formalized to keep 
the plans exercised and up-dated. It is expected that in a strategic evac¬ 
uation, the local civil defense directors will act as advisers to the normal 
government officials, rather than as commanders. Again, Hurricane Carla 
indicated that this method of operating is adequate for a strategic evacuation. 
There were no great command and control problems--even the fact that there 
was no racial segregation caused no special problems. There is a much 
greater need for planning to use available resources of trained men, working 
organizations, and materials than for the establishment of any radically new 
and different organization just to handle a strategic evacuation. 

7. 0 An Approach to Planning Assignments of Evacuees to Reception Areas 

Since Ref 3 contains a detailed description of a technique for plan¬ 
ning the assignment of evacuees to reception areas, the technique will be 
described only briefly here. The process leads first to the development of 
a suggested Distribution Plan which is a listing of the number of evacuees 
to be assigned to each reception area. Next, a Movement Plan is developed 
that tells which evacuees are assigned to each reception area, which route 
they are to use, and at what time they should depart. 
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7, 1 Distribution Plan. The distribution plan is based on the 
minimum-casualty distribution program described in Ref. 2. The fol¬ 
lowing are required items of input data: 

a. Definition of evacuation and reception sectors; 

b. Initial population, shelter capacity in each of a 
number of protection factor categories, and 
vulnerability of people in each shelter; and 

c. An Allowed Movement Table. 

The definition of evacuation areas may be made in any of a number of 

ways. However, what is basically required is a list of possible targets and 

a choice of area to be evacuated around each target. The suggested technique 

7 

makes use of the NAHICUS '63 attack study , and the area evacuated is that 
area bounded by the curve on which there is a 0. 1 probability of exceeding 
1 psi. This is the most conservative choice available if the NAHICUS results 
are used. The area evacuated would, of course, be approximated by com¬ 
monly recognized geographic features. In ' r iew of the rather large number 
of recent changes in U. S. military bases, it would b*» appropriate to make 
use of a more up-to-date attack, but the basic idea of choosing a conserva¬ 
tively large area to be evacuated can be associated with any attack picture. 

Reception areas may be chosen in a similar way. The suggested 
approach is to select areas that are aoout tne size of a county or a small 
number of counties and that have about the same vulnerability. Vulnerability 
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is defined to be the probability of exceeding some preselected casualty- 
producing radiation dose. Two hundred roentgens, the threshold for 
lethality, is suggested as a reasonable choice for this problem. 

The initial population is simply taken from census data. The shel¬ 
ter capacity is that found in the NFS Phase 2 survey. However, since the 
number of spaces in reception areas is somewhat limited, two shelter 
classes are used to represent the protection afforded by houses. A house 
with no basement is assigned a PF of 2; one with a basement is assigned a 

PF of 20. The number of basements in the reception area is estimated 

20 

from the Housing Census. Once the PF of a shelter class is chosen, the 
vulnerability of people in that shelter class is taken to be the probability of 
exceeding a free-field dose of 200 times the PF. 

The final step m preparation of input data involves the construction 
of an "Allowed Movement Table. " This table consists of an array that 
indicates whether movement between particular evacuation and reception 
areas is to be allowed or not. It is intended to help limit the problem and 
yet to allow reasonably complete use of transportation facilities. Thus, 
for example, an arbitrary distance limit might be imposed, or travel across 
a large river or other major barrier may be limited. 

20 

United States Census of Housing, i960 , Series HC(1), U. S. Department 
of Commerce, Bureau of the Census, U. S. Government Printing Office 
Janiary 1962. 
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Table IX summarizes the input data requirements for the minimum- 
casualty distribution program used to obtain the distribution plan. 

Table X summarizes the output obtained from the program. Note 
that, in addition to the distribution plan, the expected number of casualties 
for each reception sector is printed for the assumed shelter distribution 
and for the cases where the shelter system is upgraded to provide everyone 
a PF of at least 20 or 100. The number of survivors added by raising the 
minimum PF to each of these levels is also printed along with the number 
of spaces required to so upgrade the shelter system. Any other pair of 
PF's may be chosen; these were chosen to represent the value of preparing 
trenches with only weather cover (PF=20) and with about 100 psf of roof 
cover fo- shielding (PF=100). 1 

7.2 Movement Plan . Appendix C of Ref. 3 provides a descrip¬ 
tion of a computer program that may be used to design a movement plan. 

The program has since undergone considerable development, the results 

of which are reported in Ref. 4 In addition, a second technique has been 

21 

developed which allows planning of minimum-cost movements. This 

choice might be employed at a lower rung on the escalation ladder if the 

’’rate of climb" is sufficiently low. It would be appropriate, for example, to 

employ minimum-cost movements for partial evacuations of non-essential 
- ~ 

D. E. Brannon, A Computer Program fer Calculating Minimum Cost 
Movements , Dikewood Corporation Technical Note No. DC-TN-1039-6, 
December 1 7, 1964. 

For a discussion of such considerations, see Ref. 22. 
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Table DC 


Input Data for Distribution Program 

Shelter class 1 . . . n 
Capacity \'. 1 ^.nerability 

1 

2 

3 


A. Sector Characteristics 


Initial 

Sector population 


B. Allowed Movement Table 


_ To reception sectors 

From 12 3 4 


1 

2 

3 


e’ aeration 
sectors 


Table vntries :: 
movement is m 


dicate w hethen movem<■ r.t is all 
1 allowed. 1 means i* . s allo*e-' 


■ wen. ar- entrv oi 
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0 means 




Table X 


Output Data from Distribution Program 


A. Results using available she er 


Reception 

sector Shelter class 1 . . . n_ Tot al s 


1 

2 


Final population 
Expected casualties 


Final popula*'on 
Expected ^^sualties 


Total population 

Total expected casualties 


B. Results assuming minimum PF=20 or 100 


Reception 

sector 



Im PF = 20 a 


Mm PF=100 a 



Added 

Added 

Expect t d 

Added 

spaces 

Expected Added space; 

casualties 

survivors 

required 

casualties survivors require* 


2 
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Table X (Continued) 


Output Dc.ta from Distrioution Program 


C. 


Distribution table 


b 


_ Reception s ctor 

Eva cuation 

sector 12 3 


1 

2 


3 . 

Any pair of PF's may be used. 

Table entries are numbers of people moved from an evacuation sector 
to a reception sector. 
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workers and their dependents, particularly if the crisis were still at levels at 
which decision-makers would not want to disrupt the transportation system. The 
input and output associated with each of these programs is described below. 

7.2.1 A Computer Program for Planning Rapid Strategic Movement . 
The compuier program described in Ref. 4 operates on an initial population 
distribution, a desired final distribution, movement rates and a segment-node 
description of the route network over which the movement is to be accom¬ 
plished. From this data, the number of people to be moved from each evac¬ 
uation site to each reception site, the routes involved, and the associated 
time schedules are calculated. The logic in the program is basically heuristic 
and consists of a series of algorithms, originated during the development of 
the program, that tend to minimize the time required to attain the desired 
distribution. 

A mathematical proof that the technique used leads to a minimum¬ 
time movement has not been found. Attempts to develop a technique that 
has a more rigorous basis have been unsuccessful to date. 

The initial population is taken from census data and should corre¬ 
spond to that used in the preparation of input for the Distribution Program 
(see Table IX). The desired final distribution is taken from the output of 
the Distribution Program (Part C of Table X), Note that the Distribution 
Program output associates a number of evacuees with each evacuation- 
reception sector pair. However, it is the desired final distribution that is 
of greatest interest, since any other association of evacuees with evacuation- 
reception sector pairs tnat gives the same final distribution will also be a 
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minimum-casualty movement. The Movement Program is therefore em¬ 
ployed to provide a set of such associations that will make efficient use of 
the transportation system. The present version of this Movement Program 
will accept as many as 10 evacuation sectors and 50 reception sectors. The 
Distribution Program will treat many more, but if the Movement Program 
limits are exceeded, the Distribution Program output will have to be broken 
into parts for input to the Movement Program. 

The evacuation rate is defined in terms of the number of people per 
hour that can move past a point on a "unit-capacity" segment of a trans¬ 
portation route. A unit-capacity segment may, for example, be defined 
as one lane of a highway. Then, if one assumes an evacuation rate of 1000 
automobiles per lane per hour, with each auto containing 4 persons, the 
evacuation rate is 4000 people per hour. Any other physically reasonable 
choice can be used. 

A general speed in miles per hour is specitieo that is consistent 
with the probable travel rate over the evacuation routes (e. g. , 30 mph). 

If there is a segment over which a speed can be maintained that is signif¬ 
icantly different from the general speed, the mileage of this segment must 
be appropriately adjusted. Thus, if 4 5 mph can be maintained on a given 
segment while the general route speed is 30 mph, the true segment mileage 
is multiplied by 30/45 or 0. 67 to adjust its length. This kind of adjustment 
is also used to treat othe 1 " types of transportation, such as railroads. Here 
the adjustment can be made to reflect necessary changes in both evacuation 
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rate and route speed, the first by varying the number of unit-capacity seg¬ 
ments, and the second by adjustment of segment length. 

Memory of the Dikewood 7044 computer limits the present version 
of the program to treatment of less than 500 segments. A few remarks to 
explain the node-numbering conventions used in the program will help the 
reader understand the sample problem discussed later in this section. How¬ 
ever, these remarks are not essential to understanding the program and its 
limitations and the reader not interested in these details should skip to the 
text that follows Fig. 6. The node numbers assigned to evacuation areas 
must be less than 100, and the numuers assigned serially to reception areas 
must start at 2000. Non-terminal nodes are numbered serially beginning with 
100 and numbers from 1300 to 1399 are used for dummy nodes to represent 
multiple-capacity route segments. Some reception places will also be junc¬ 
tions in the network for travel to other reception sectors. Such places are 
assigned a non-terminal node number for the travel network and a zero- 
length segment with a terminal node number to indicate that it is a recep¬ 
tion site. 

Figure 4 shows a single 10- mile segment of a route between nodes 
220 and 221. 



Map of a Unit-Capacity Route Segment 
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If a portion of a route has a capacity of two (two lanes of a highway 
to be used) it is mapped with a dummy node into two unit-capacity segments. 
Figure 5 shows a 12-mile, 2-lane segment available between nodes 221 and 
222 . 



Fig. 5 

Map of a Double-Capacity Segment 

When a route to a given reception area passes through another recep 
tion area, a dummy node and a zero-length segment are used so that the 
bypassed reception area remains a terminal node. Figure 6 illustrates this 
situation. 



Map of Route to Reception Area 2002 that Passes 
Through Reception Area 2023 
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The direction of traffic flow must also be indicated in the input data. 

In Figs. 4, 5, and 6, this direction is indicated by the placement of an arrow¬ 
head. This choice requires considerable judgment on the pari of the planner. 
He should first make some gross estimates of the desired overall direction 
of flow. In making the final choices, he will find that a large problem area 
may be divided into several subareas simply by proper choice of flow direc¬ 
tions on individual segments. 

Two additional input data items must be specified, namely, a time 
unit and a precision variable. The time unit is specified for the problem 
such that the physical length of the line of people passing a point in unit 
time is small compared to the number of people using the route. For 
example, a time unit of 0. 01 hours and an evacuation rate of 4000 people 
per hour would establish the program's concern with the movement of units 
containing 40 people. The precision variable establishes the degree of 
precision desired in the calculations and is equal to the desired precision 
divided by the time unit. For example, if the desired precision is a half- 
hour (that is, no readjustments are to be made in the schedule if it cannot 
be improved by more than a half-hour) and the time unit is 0. 01, the pre- 
cision variable is 0.5/0. 01 or 50. 

The output of the program is simply a list of the numbers of people 
to be moved from each evacuation site to each recepUon site, the route to 
be taken by each group, the time of departure after movement begins, ar.d 
the time when each reception sector is full. 




Several sample problems were run in the development of the pro¬ 
gram, Among these was one used to produce a sample strategic movement 

plan for the State of New Mexico. Target cities within the state and in 

g 

neighboring states were taken from the Tech-Ops Target List. It was 
assumed that half of the population of El Paso, Texas, would also take part 
in the New Mexico movement. This was done because of the scarcity of 
reception areas, external to New Mexico, for the El Paso population. All 
other places within the state with a population greater than 200 people were 
divided into groups to form reception areas Each reception area was given 
the name of the largest place within it. 

Figure 7 is an arrow diagram map of allowed traffic flow super¬ 
imposed on a segment-node map of the highway network. 

The number of people to be sent to each reception area was based 
on a uniform load factor for the state. (The load factor was 2. 28 for the 
problem. ) Table XI lists the initial population of each reception area and 
the number of evacuees assigned to it. 

Figure 8 depicts the evacuation routes that were calculated by the 
program. Evacuation areas are enclosed by 24-mile radius circles. 
Principal cities in reception areas are designated by crosses. The com¬ 
puter printout of the schedule is shown in Table XII. Note that there are 
6 evacuation areas and 38 reception areas. There are 42 different routes 
used in the schedule and 27 hours are required to complete the movement 
of the 637, 616 people by auto. No rail movement r > were used in this problem. 
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Segment-Node Diagram for New Mexico Strategic Movement 



Table XI 
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Evacuee Assignment for a New Mexico Strategic Movement Plan 


Reception area 1 

Initial 

a 

Dopulation 

Number of 
evacuees assigned 

Farmington 

53, 306 

68, 348 

Gallup 

37,209 

47, 709 

Grants 

22,939 

29,412 

Cuba 

5,469 

7, 012 

Tierra Amarilla 

7, 443 

9, 543 

Vallecitos 

958 

1, 228 

Taos 

24, 653 

31, 610 

Raton 

10, 408 

13, 345 

Clayton 

6, 068 

7, 780 

Springer 

3, 398 

4, 357 

Mosquero 

1,875 

2, 404 

Logan 

674 

864 

Turumcari 

11,605 

14, 880 

Las Vegas 

23, 468 

30, 090 

Santa Fe 

38,388 

49, 220 

Santa Rosa 

4, 308 

5, 524 

Moriarty 

3, 073 

3, 940 

Belen 

9, 101 

11, 669 

Mountainair 

2, 356 

3, 034 

Encino 

1, 058 

1, 356 

Vaughn 

1, 302 

\ 669 

Jal 

8 927 

11, 446 

Reserve 

2, 773 

3, 555 

Magdalena 

1, 825 

2, 340 

Socorro 

8, 343 

10, 697 

Ft. Sumner 

2, 991 

3, 835 

Carrizozo 

2, 571 

2, 296 

Hondo 

5, 173 

6, 633 

Truth or Consequences 

6,409 

8, 217 

Silver City 

21,059 

27, 001 

Dora 

2, 781 

3, 566 

Artesia 

1 7, 686 

22,677 

Lovington 

15, 034 

19, 276 

Carlsbad 

33, 397 

42,821 

Hobbs 

29, 468 

37, 783 

Lordsburg 

5,215 

6, 687 

Deming 

9, 839 

12, 615 

Las Cruces 

54,728 

70 ’77 

Totals 

497, 288 

637, 61 6 

£ 

Population based on 1960 census 



Includes 205, 350 evacuees from 

El Paso. 
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Fig. 8 


of a Strategic Movement Plan for New Mexico 
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ROUTE NUMBER 113 ALAMOGORDO TO SOCORRO - MOVE 10413 PEOPLE WITH START TIME 

ALAMOGORDO TO CARR120Z0 TO SAN ANTONIO TO SOCORRO 
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Another problem for which this Drogram was used is described in 
Par. 2. 2 of Chapter III. 

7.2.2 A Computer Program for Planning Minimum- Co st Movements. 
A second movement planning technique has been prepared for possible use 
following strategic warning. Emphasis is placed on the word strategic; it 
is assumed that the decision-maker has high confidence that the warning is 
strategic and that he would therefore prefer to carry out a minimum-cost 
movement. Cost is expressed in terms of man-miles and the program 
determines that movement which minimizes total man-miles traveled. It 
is further assumed that only one lane of each route would be made avail¬ 
able for evacuees. This is considered appropriate for the assumed level 
of crisis; a level at which one may, for example, wish to evacuate depend¬ 
ents, but not preclude use of the transportation system for other essential 
purposes. Keeping other lanes open would also simplify logistics and con¬ 
trol requirements. 

The input data required for this program is very similar to that for 
the program described in Par. 7. 2.1. A principal difference is that this 
program does not require prior choice of direction of flow and the prepa¬ 
ration of the node-segment network is therefore much simpler. Otherwise, 
the program requires initial and final population distributions, evacuation 
rate, and the general speed of travel ever the transportation network to 
be specified. 
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The calculations are also more straightforward, First an algorithm 
23 

by Moore is applied to calculate the shortest route from each evacuation 
site to each reception site. These shortest routes are stored in the form of 
a sequence of node numbers and total route length. The array of shortest - 
route lengths between each evacuation site and each reception site may be 
thought of as a "cost" matrix. 

24 

The next step in the calculations makes use of another algorithm 
and the initial and final population distributions and the cost matrix. This 
algorithm allows calculation of the number of people to be moved from each 
evacuation site to each reception site such that the total number of man- 
miles traveled is a minimum. 

In the next section of the program, competition is resolved on those 
routes selected for use. Evacuees moving the greatest distance are given 
first priority. If two groups of evacuees both use a particular segment of 
highway, the group traveling the shorter distance is delayed at their ori¬ 
gin long enough so that the segment in common use is available to them 
just as they reach it. 

23 

Pollack, Maurice and Wiehenson, Walter, "Solutions of the Shortest- 
Route Problem--A Review, " Operations Research , Vol. 12, No. 4, 
pp. 519-653. 

24 

Ford, L. R., and Fulkerson, D. R., Flows in Networks , Princeton 
University Press, 1962, p. 95. 
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The other choice that would use the highway efficiently would require 
the group traveling the shorter distance to wait in a "holding area" near the 
entrance to the segment in common use. Such a procedure may be worth¬ 
while in movements carried out at a higher level of crisis to get people out 
of potential initial effects areas as soon as possible. However, for lower- 
level crises appropriate to choice of a minimum-cost movement, asking 
evacuees to wait in holding areas for times as long as a day seems undesirable. 

Finally, the program prints a listing of the route description, 
departure schedule, and total time required to complete movement to each 
reception site. 

This program has also been applied to a number of problems, oi*c 
of which is described here to illustrate use of the program. The place to 
be evacuated is Albuquerque, New’ Mexico; the reception places are arbi¬ 
trarily chosen to include all places with populations over 200 in about the 
northern half of the state. The input data assumptions are listed below. 

The initial distribution of population was taken from the 1960 Census. 

The desired final distribution of evacuees was found using the minimum- 

casualty movement technique developed by Dikcwood under Contract 

1 2 

OCD-OS-62-248. ' As indicated previously, this technique requires 
assumptions concerning the attack, shelter distribution, and a casualty 
criterion. The NAHICUS ’63 results were used to describe the fallout 
threat in reception areas the the area evacuated was a 24-mile radius 
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circle around the aiming points It was assumed that NFSS Phase I shelters 
and mines from the Phase II survey were used for fallout protection A 
PF of two, associated with houses, was assigned to those for whom there 
is no space in the NFSS shelters. Mortality was selected as the v.dsualty 
level to be minimized; a simple straight-line mortality-versus-dose rela¬ 
tion was assumed, with the threshold at 200 roentgens and 100 percent 
mortality at 750 roentgens. A maximum housing load factor of three was 
chosen and a calculation made of the desired distribution of evacuees. The 
results are shown in Table XIII. 

The desired distribution shown in Table XIII was then used as part 
cf the input date! to the transportation assignment program described here. 
The other input required for this calculation is a node-link diagram of the 
transportation system. This diagram was prepared for the Albuquerque 
area assuming that one lane of each mailable highway would be used; the 
results of the transportation ac,i i j...ncr't program are given Table XIV. 
Figure 9 consists of a map of the highway network used as input to the 
program. Fig 10 illustrates the solution obtained. 

Table XV consists of a printout of the computer output for the sam¬ 
ple problem 
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Table XIII 


Desired Population Distribution 
for Minimum-Fatality Movement from Albuquerque 
(Maximum load factor * 3) 

(a) 

Population (thousands) 


Sector No. 

Principal city 

Initial 

Final 

1-6 

Albuquerque 

274. 6 

0 

7 

Farmington 

31. 1 

93.0 

8 

Santa Fe 

37. 6 

108. 7 

9 

Gallup 

14. 9 

14. 9 

10 

Grants 

14. 7 

44. 1 

11 

Belen 

5. 6 

16. 8 

12 

Clayton 

4. 5 

13. 5 

13 

Las Vegas 

15. 1 

45. 3 

14 

Raton 

8.4 

25. 1 

15 

Socorro 

6.2 

18. 6 

16 

Magdalena 

1.4 

4.2 

17 

Mountamair 

2. ' 

8. 1 

18 

Taos 

6. 7 

20. 1 

19 

Vaughn 

1. 5 

1.5 

20 

Santa Rosa 

2. 5 

2. 5 

21 

Springer 

4. 0 

12.0 

22 

Cuba 

1. 5 

4. 6 


a 1960 Census. 
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CHAPTER III 


EVALUATION OF STRATEGIC MOVEMENT PLANS 

1.0 Introduction 

Some techniques for planning strategic movements were discussed 

in Chapter II. These techniques would produce plans based on movements 

that minimize expected casualties in the face of a heavy attack. This chap- 

* 

ter is concerned with the evaluation of such planned movements under 
various particular situations. For example, one may wish to evaluate the 
effects of some particular wind structures, of various attacks, or of vari¬ 
ous alternative responses when an attack arrives before a movement is 
completed. Standard target-analysis techniques can be applied to the first 
two of these questions, but the third required development of a new tech¬ 
nique. Again, the computer was found helpful because the bookkeeping 

problem gets quite involved. The computer program developed for this 
25 

purpose requires as input, the number of survivors among those left in 
the target area at the time of attack. This calculation of survivors is 
easy to do by hand, but the volume of work foreseen seemed to indicate 


*The evaluation is expressed in terms of survivors versus time to attack, 
as in curves A and B in Fig. 1. 

25 

D. E Brannon, A Computer Program for Calculating Fatalities Among 
Evacuees Enroute When an Attack Begins , Dike wood Technical Note No. 
DC-TN-1039-5, December 3, 1964. 
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that it . 'uid be useful to program this problem also. This chapter con¬ 
sists primarily of a description of these two programs. A few applications 
of the programs are also described. 


2. 0 Casualty Calculations for Evacuees Enroute When an Attack Begins 

2.1 Description of the Computer Program. As indicated in the 
introduction to this chapter, a computer program has been developed that 
allows calculation of casualties among those enroute as well as among 
those in either a target or reception area at the time of attack. The pro¬ 
gram accounts for changes in location and protection factor after the attack 
as well as before. 

The following input data items are required for this program: 

1. Initial and desired final distribution of evacuees. 

2. Number of evacuees assigned to each route. 

3. Schedule of evacuee departures, rates of movement, 
and protection factors. One intermediate stopping 
point is allowed. 

4. Latitudes and longitudes of enough points to de¬ 
scribe each route. 

5. Weapon yields, fission fractions, burst locations 
(latitude, longitude, and surface or airburst), and 
detonation times after movement begins. 

26 

D. E. Brannon, A Computer Program for Calculating Casualties From 
the Initial Effects of Nuclear Weapons , Dikewood Technical Note No. 
DC-TN-1039-8, February 4, 1965. 
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6. Effective wind velocity and shear components. 


27 

7. Fraction of those in the target area who become 
casualties from initial effects. 

Census data is usually a satisfactory source for the initial population 
distribution. The desired final distribution, route assignments, and move¬ 
ment schedules may be obtained using the procedures outlined in Chapter II. 
The protection factor and time unsheltered after arrival in the reception 
area must be specified; only one value of each of these two variables is used 
to represent all reception areas. One intermediate stopping point is allowed 
on each route. This feature was incorporated to permit analysis of the value 
of "holding areas. " Holding areas can be large parking lots located, for 
example, on the city side of the point where the road leaving a target city 
narrows down from four to two lanes. Latitudes and longitudes oi points 
along a route (Item 4) may be obtained from a number of types of maps; 
the distance between points need not be kept constant. The attack assump¬ 
tions required (Item 5) need no further explanation. Determination of the 
effective wind speed and shear components (Item 6) is fairly complicated, 
but is described in Ref. 27 which also contains a description of the fallout 
model used as a subroutine in this program. A computer program has been 
prepared to convert observed wind structures to the desired form. Anyone 

27 

Wood, W. D., et al.. Emergency Operations Doctrine and Organization . 

Dike vood Corporation Report No. QR-1040-2, Addendum No. 1; 

February 14, 1964. (Secret) 
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who has many problems to solve may find it useful; copies may be obtained 
by writing to one of the authors of this report. As indicated in Par. 1. 0, 
the fraction of thooc in me target area who become casualties from initial 
effects (Item 7) may be calculated using standard target-analysis techniaues 
or the computer program described in Par. 3. 0 of this chapter. 

To calculate casualties, the computer is used to simulate the move¬ 
ment of people along each route. The events that take place and their asso¬ 
ciated time periods are illustrated in Fig. 11 and listed below. 


t=0 


a 



Fig. 11 


t=2 weeks 
after last 
detonation 


History of Events for Enroute Casualty Calculations 

a. The delay time before using the route (experienced by all 
units using the route). 

b. Further delay experienced by a unit while all of the preceding 
units start out over the route. 

c. The time period required for any unit to travel to a point of 
intermediate delay or "holding" area. 

d. The delay experienced in the holding area. 

e. Travel time from holding area to reception area. 

f. The time period for which units are unsheltered after arriving 
in the reception area. 

g. The time period for which each unit is sheltered. 

T. Time of onset of attack (may occur at any assumed time). 
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During time periods (a) and (b) the users of a route may accumulate 
dose from fallout and may become fatalities from initial effects, if any are 
experienced at the coordinate position denoting the starting point for the 
route. 

Figure 12 indicates the layout of a route, and shows the method of 
dose calculation for a unit of people. 


Unit of people 

czrz: i: 


A 


'✓/✓// 



Column of people 
using route 


$-f-©H-1—- \e h r » 

A N-lI N 

I 

P 



♦ o j— i—4> 

B 

o --Route descriptor 
points 


Fig. 12 

Illustration of Technique Used for Dose Calculations 

For example, consider the calculation of the dose added when the 
leading member of a unit advances from the point (N-l) to N, where N is 
an integral number of distance increments (DELL) from the origin. The 
dose rate at the coordinates midway between N-l and N (point P) and the 
time t is calculated. The additional dose received by the unit at P is 
then given by: 

A dose * dose rate (P, t)At, 
where At is the time unit (DELL/movement speed). 
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If the movement plan includes a stop at a holding area, the dose is 
accumulated at the point chosen to represent the holding area for the delay 
period represented by (d) in Fig. 11. 

For times after arrival in the reception area, the dose received by 
a unit is calculated for the fixed coordinates of the reception center. After 
period (f) has passed for a unit, the people are assumed to be sheltered with 
the specified protection factor out to a time of two weeks after the last weap¬ 
on detonation. The two-week dose is considered a reasonable maximum, 
because movement, decontamination measures, and radioactive decay will 
help reduce the exposure to relatively low levels after two weeks have passed. 
However, if desired, it is a simple matter to change the program to make 
use of any other upper limit of integration. 

Direct numerical integration of the dose rate to find total dose is 
performed for a period of 24 hours after the last burst, to properly account 

for the dose accumulated during fallout arrival. From then on, the fallout 

- 1.2 

is assumed to be on the ground and decaying according to the t law. 

The effect of biological recovery is ignored in the present program. 

The number of fatalities from fallout experienced within a unit of 
people is calculated from a straight-line approximation of deaths versus 
dose between the limits of 200 r (0% fatalities) and 750 r (100% fatalities). 
Further, if thr happens to be within a 24-mile radius of a target point 
when the target is attacked, the fatalities in the unit from primary effects 
are also calculated as a fraction of the number of people in the unit. 
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The fatalities for all units are added together to find total fatalities for those 
using a given route. 

The printed output from the program includes the items indicated in 
Table XVI. 

Table XVI 

Output Information for Enroute Casualty Calculation Program 

1. Attack and shelter variables. 

2. For each route: 

a. Origin and destination, 

b. Route number, 

c. Planned number of people to be moved over the route, 

d. Number of casualties, 

e. Percent casualties, and 

f. Dose xeceived by each unit of people using the route. 

3. For each problem: 

a. Planned total number of people to be moved over all 
routes, 

b. Total casualties, 

c. Total percent casualties, 

d. People per unit, and 

e. Time from beginning of movement to first detonation. 

2.2 A Sample Problem . As indicated in Par. 2.1, this program 
requires as input, a schedule of evacuee departures and rates of movement. 
To obtain this schedule, the problem used in Par. 7.2. 2 of Chapter n to 
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illustrate the minimum man-mile movement program was run using the 
other movement program described in Chapter II, Par. 7. 2.1. The initial 
and desired final population distributions are given in Table XIII in Chap¬ 
ter II. The movement schedule obtained is listed in Table XVII and is 
depicted in Fig. 13. The PF was assumed to be unity until one hour after 
arrival in the reception area at which time PF a 2. This is equivalent to 
assigning evacuees the protection factor afforded by a house. 

The following attack conditions were assumed: 

Yield * 10 Megatons 

Fission fraction * 0. 5 
Burst height * surface 

Burst location * intersection of major runways 

at Kirtland Air Force Base 

Burst time * 8 hours after movement begins. 

A wind structure observed at Albuquerque on 12 May 1956 was used; 

the values used xn the computation were: 

Effective windspeed * 50 knots 

Effective wind angle ■ 45 degrees 

Downwind (x) shear component ■ ♦O. 15 knots/kft 

Crosswind (y) shear component « -0.027 knots/kft 

A hand calculation was made which indicated that about 60 percent 

of the population remaining in Albuquerque at the time of attack would be 

killed from inrial effects if they were distributed as in the 1960 Census. 
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Notes: 

1. Principal places in reception areas are 
shown with no. evacuees thousands), 
least start time delay (hours), and 
time when full, respectively. 

2. New Mexico highway numbers are in 
circles; US in rectangles. 



Fig. 13 
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For illustrative purposes, the surviving evacuees were assumed to 
move according to the preplanned schedule after the attack takes place. 

The detailed results are shown in Table XVIII. 

Among other things, the results in Table XVIII indicate that the assumed 
conditions lead to a total of 37 percent fatalities among the Albuquerque popu¬ 
lation if the attack occurs eight hours after movement begins. As indicated 
in Table XVI, the list of numbers following the results for each route show 
the two-week dose received by each unit of evacuees. For example, if those 
going to Santa Fe are to survive, protective measures that would reduce their 
doses by an additional factor of about six are required. 

3. 0 Calculation of Casualties from Initial Effects 

3.1 Description of the Initial Effects Casualty Calculation Program. 
This program is used to calculate casualties within a target city for various 
attack and shielding assumptions. The calculations are made under the 
assumption that the vulnerability of the target city population can be repre¬ 
sented by the following shielding categories, the first seven of which are 
defined in Ref. 28: 

1. Outside unshielded 

2. Outside shielded 

3. Wood frame structure 

2 3 

Davis, L. W., et al., Pr ediction of Urban Casualties From the Immediate 
Effects of a Nuclear Attack , Contract OCD-OS-62-203, Dikewood Report 
No. DC-FR-1028; April, 1963. (Confidential) 
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ROUTE NUMBER PEOPLE MOTED EXPECTED CASUALTIES PERCENT CASUALTIES 
142 12360.00 0. 0. 



ROUTE NUMBER PEOPLE MOVED EXPECTED CASUALTIES PERCENT CASUALTIES 
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4. Brick structure 

5. Light steel frame structure 

6. Heavy steel frame structure 

7. Reinforced concrete structure 

8. Undefined 

9. Undefined 

10. Undefined. 

The last three categories were left open and undefined to allow for 

other shielding postures that one might make use of in a particular problem; 

e.g., for people in automobiles, basements, trenches, or blast shelters. 

The mortality curves associated with each shielding category have been 

converted to tabulated values of the fraction of people killed as a function of 

distance, burst height, and yield. Tables for two values of burst height 

art available for the first 7 categories, namely, zero and the Hiroshima 

1/3 

scaled height of burst, 846W . The Hiroshima height is optimum for 

29 

about 8 psi. Tables are available for yields of 0. 4, 1, 4, 10, 25, and 
50 Mt. 

The following input data is also required: 

1. Number of weapons, their yield, burst 
height (surface or air), and burst point 
latitude and longitude. 

29 • 

Samuel Glasstone, Editor, The Effects of Nuclear Weapons , Supt. of 
Documents, U. S. Government Printing Office, Washington, D. C.; 

1962, ($3.00) 


- 78 - 


2. Geographic unit number, e. g., Standard 

30 

Location Area (SLA) number, and a sin¬ 
gle latitude and longitude pair used to repre¬ 
sent the location of the people in that unit. 

3. Number of people in each shielding cate¬ 
gory in each geographic unit. 

The program simply calculates the distance from the burst point to 
each geographic unit, look? up the fraction of survivors in each shielding 
category, and multiplies by the number in that category. For multiple 
weapon attacks it is assumed that the weapons act independently, i. e., the 
survivors of the first explosion are assumed to be in the same shielding 
category after the first explosion as before. 

The program output includes: 

1. City name 

2. Weapon yield, burst height, and burst 
location 

3. For each geographic unit: 

a. Latitude and longitude 

b. Peak overpressure 

c. Survivors in each of the ten shielding 
categories 

d. Percent survivors m the geographic unit 

4. Total population, total number, and total 
percent aurv^-ors 

30 

National Location Code , OCD-OEP Region 5, OCD-FG-IV-3. 115, Bureau 
of the Census, 1962. 
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3. 2 A Sample Initial Effects Casualty Calculation. The attack 


assumptions for this sample problem are the same as for the sample prob¬ 
lem in Par. 2.2 of this chapter. 

Standard Location Areas »vere taken as the basic geographic units. 
The population in each SLA was taken from Ref. 30. To obtain the distri¬ 
bution of the people in each SLA over the ten shielding categories, estimates 

were made of the type of building the people were housed in by examining the 

20 

United States Census of Housing for 1960, and Sanborn Maps of the area. 
These estimates indicated that only 3 of the 10 shielding categories were 
occupied assuming that the population was unwarned. These three categories 
were: wood-frame structures, brick structures, and reinforced concrete 
structures (apartment houses). 

The results of the calculation are reproduced in Table XIX. Note 
that 56 percent of the people are killed, which is in reasonable agreement 
with the hand-calculated value of 60 percent used in Par. 2.2. The hand- 
calculated value is somewhat less accurate because it was based on a less 
detailed examination of the shielding posture of the population. 

4. 0 Some Implications of Various Trans-Attack Policies for Strategic 
Movements Interrupted by War 

A number of runs of the program described in Par. 2.0 have been 
made to estimate casualties for a variety of input parameters and trans¬ 
attack policies for strategic movements interrupted by war. 
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At the time the calculations discussed in this paragraph were made, 
the Movement program (Chapter II, Par. 7. 2. 1) was nearing completion 
of its development. The available version of the program was applied to 
the problem described in Par. 2. 2 of this chapter and the results obtained 
were not greatly different from those given in Table XVII and Fig. 13. 

Thus, while the results to be described at this point are based on a slightly 
different movement, the reader can get a sufficiently accurate picture of 
the movement schedule from Fig. 13 

The assumed attack consists of a 10-Mt, 50-percent fission weapon 
surface burst at a runway intersection at Khrtlanu Air Force Base in 
Albuquerque. 

The wind velocity and shear components were determined from wind 

g 

structures used to develop Fig. 2. These structures are believed to repre¬ 
sent, reasonably well, the range of structures to be expected in the Albu¬ 
querque area. 

It was assumed that the evacuees would be unshelteied for one hour 
after arrival in their reception place, after which they would be given a 
PF of two. It was assumed that the movement would continue as planned 
after the attack took place. 

Figure 14 illustrates the results lor the four wind structures. It 
shows that more than half the evacuees can be expected to su^uve if the 
attack takes place as early as seven hours after the movement begins. If 
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Surviving fraction of evacuees 




the attack occurs just as the movement is completed, a minimum of 0. 7 of 


the evacuees would survive compared to a maximum of 0. 3 if the attack 
occurs just as the movement begins. 

Figure 14 also illustrates the effect of using holding areas near the 
edge of the initial effects area. For this one set of calculations, an improve¬ 
ment of about 15 percent in the fraction surviving is realized for attacks 
tha, occur about 4 hours after the movement begins. 

Figure 15 shows the effect of employing simple trench or basement 
shelters. Note that, if the movement is completed, all evacuees would sur¬ 
vive and even if the attack occurs as the movement begins, there is about a 
30 percent increase in the fraction surviving. Of course, this is optimistic 
because of the assumption that, after the attack, the survivors would continue 
the movement as planned. 

Figure 16 shows the effect of a 24-hour delay in the construction of 
shelters; 24 hours is a conservatively large estimate of the time required 
to prepare a trench snelter. Since the gain is small if the evacuees are 
unsheltered 24 hours, it would be worth-while to begin shelter construction 
at a lower rung of the escalation ladder than the one used to Degin the final 
stage of an evacuation. 

Some more detailed calculations along these lines seem in order. 

Of special interest would be the calculation of casualties when the population 
location and vulnerability in town are improved by moving people out of the 
area3 near the likely aiming points before moving others. This would re¬ 
quire development of a detailed movement plan within the city, but it would 
appear to be worth the effort. 
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Time from beginning of movement to attack (hours) 

Fig 15 

Survivors Added by Improving Protection Factor 
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Surviving fraction of evacuees 
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CHAPTER IV 


SOME SHELTER POLICIES AND PROGRAMS 


1.0 Introduction 

This chapter consists of a discussion of the development of curves of 
Type C in Fig. 1. Once again, as for movement strategies, one may wish to 
first plan a shelter system, then to evaluate it against a range of attacks. 

The evaluation against a particular attack to obtain curves of Type C is a 
relatively simple matter, but the shelter system design may not be. A num¬ 
ber of shelter system design policies have been suggested by various authors. 

The following are among the policies considered of possible interest: 

* 

1. Uniform maximum fatalities per incoming weapon. 

2. Uniform probability of survival. 

3. Uniform shelter design overpressure. 

4. Maximum added survivors per dollar expended. 


5. Maximum enemy cost per kill. 

The first policy was originally suggested by Hudson Institute (Herman 

31 

Kahn and William Brown). Basically, the policy leads to a shelter system 
that gives up a constant maximum number, 0, of fatalities for each incoming 
weapon. Shelter hardness is then related to population density, anticipated 
yields of the incoming weapons, budget, and 0. Table XX shows how these 


Note added in proof: A recent IDA report, S-186, describes a computer 
program for Policy (1). 


W. M. Brown, The Design and Performance of "Optimum" Blast Shelter 
Programs , Hudson Institute Report No. HI-361-RR/2; June 11, 1964. 
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Table XX 


Cost of Blast Shelters for 213 Urbanized Areas (1960) 


Criterion, 0 
(thousands 
of fatalities) 

Shelter required (psi) at 

Total 

cost 

(billioi 

p =3, 000 

p = 10, 000 

p = 30, 000 

5 

90 

300 

900 

31.6 

10 

45 

150 

450 

23. 8 

15 

30 

100 

300 

20. 3 

25 

18 

60 

180 

16. 8 

40 

11 

37 

112 

14. 3 

60 

7. 5 

25 

75 

12. 5 

80 

5. 6 

19 

56 

11. 5 


❖ 

p = population density (number per square mile). 
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variables are related for some sample values. The table was constructed 

associating one-megaton weapons with the mortality criterion, 0. 

This policy seems appropriate for a shelter system designed against 

a population attack. However, other policies may also prove interesting for 

various reasons. For example, if a counterforce target is located in the 

vicinity of a city, it may be appropriate to consider the influence of the target 

on the choice of shelter hardness for that city. Also, it seems appropriate 

to compare systems based on the various policies to determine whether some 

conclusions can be drawn that are insensitive to the choice of policy. Or, 

Policy 4 in the above list might be used to obtain a preferred order of 

* 

spending money on whichever ot'.er policy is chosen as a design basis. 

Computer programs for designing systems associated with Policies 
2, 3, and 4 have been prepared and are briefly discussed in this chapter. 

If the reader is interested in using the programs, a detailed description is 
provided in Refs. 33. 34, and 35. 


32 

Herman Kahn, Some Comments on Group A, B, and C's Work. A Project 
Harbor Briefing . Hudson Institute Report No. HI-305-BN, November 15, 1963. 

lit 

The computer program associated with Policy (4) can also be used to obtain 
a preferred order of spending for other protective programs, such as thermal 
countermeasures. 

33 

K. D. Granz jw and D. L. Summers, A Computer Program for Calculating 
Shelter Hardness and Cost for Uniform Survival Probability. Dikewood 
Corporation Technical Note No. DC-TN-1039-3; December 3, 1964. 

34 

K. D. Granzow, et al., A Computer Program for Calculating the Expected 
Cost per Survivor for a Uniform Shelter Design Overpressure Policy . 
Dikewood Corporation Technical Note No. DC-TN-1039-7; January 29, 1965. 

35 

D. L. Summers and A. R. Bliss. A Computer Program for Finding the Order 
of Expenditure on Blast Shelters that Maximizes Additional Survivors per 
Dollar . Dikewood Corporation Technical Note No. DC-TN-1039-10; June 1965. 
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The expected number of fatalities could be calculated for use in 
Policy 5 studies with a simple modification of a subroutine used in the 
three programs associated with Policies 2, 3, and 4. This has not been 
done, however, and the rest of this chapter is concerned only with an 
explanation of the computer programs for Policies 2, 3, and 4. 

Policies 2 and 3 are fairly obvious possibilities and had been se¬ 
lected for study when it was learned that W. F. Roherty at Sandia Cor¬ 
poration had already developed computer programs to calculate the impli¬ 
cations of Policy 4 in addition to 2 and 3. He discussed his work with 
Dikewood staff members and, while the programs have undergone several 
revisions, they still make considerable use of ideas and actual routines 
originally suggested by Mr. Roherty. The principal changes include a sub¬ 
routine to calculate either the geometric or the population centroid of an 
area of interest and the development of a table look-up procedure for 
finding the probability of covering a point a known distance from an aiming 
point with at least a given initial effect. The table look-up procedure is 
faster than the calculation procedure used in the original program. Provi¬ 
sions have also been made to include cost functions other than the one in the 
original program. Other changes are principally in format and similar 
details. 

One further acknowledgement is in order, namely, to Mr. Luke 
Vortman also of Sandia Corporation. Mr. Vortman had applied the original 
programs tc perform an analysis of shelter systems for Albuquerque and he 
made his results available to Dikewood for use in checking out the new pro¬ 
grams. 
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2. 0 Uniform Probability of Survival 


This program was written to facilitate calculation of the shelter design 
overpressures required to give a uniform survival probability to the population 
considered. The area of interest is divided into tracts small enough to permit 
the effects of assumed weapons to be considered uniform over the tract. Al¬ 
though SLA's do not always satisfy this condition, their use is believed appro¬ 
priate since uncertainties about the enemy choice of aiming point and yield 
probably override the errors associated with using SLA's. The position 
(lat. and long. ) and the population of each tract are input data to the program. 
A table of survival percentages (up to ten at a time) for which the calculation 
is to be done must also be supplied. 

The program is arranged so that a specific aiming point can be given 
as input or the program will compute the geometric or population centroid of 
the tracts and use either one of these as the aiming point. A CEP and a table 
of radius of effect versus shelter design overpressure are also required 

Besides yielding a shelter design overpressure for each tract for each 
survival percentage, the program computes the total cost and the cost per 
survivor for the shelters associated with each survival percentage for the 
overall area of interest. These costs are computed in three different ways 

¥ M 

(1) by the formula cost ($)per person* 3(Wdesign overpressure(psi) 

Originally suggested by Mr, Luke Vortman, Sandia Corporation, but 
obviously not intended for use at low overpressures 


-92- 



3l 

(2) by the formula: cost ($) per person * 50+2 (W~de sign over pres sure (psi) 

(3) by using a table of cost versus design overpressure supplied as 
input to the computer. 

A table of values is used to relate the probability of covering a point a 
given distance, 6 , from the aiming point with a given overpressure (radius of 
effect). The tabic, lists probability of coverage as a function of radius of effect 
divided by CEP and 6 divided by CEP. 

For a particular tract or SLA, the program finds the probability of 
coverage with each of the overpressures which are listed in order of in¬ 
creasing overpressure. A percent-saved figure is then chosen, converted 
to a decimal and subtracted from one to give the corresponding mortality 
fraction. This fraction is then compared to each of the probability of cover¬ 
age figures to find the first overpressure at which the probability of coverage 
is less than or equal to the value of the fraction. This implies that at this 
overpressure at least the desired percent of the population is saved, i. e., 
the mortality fraction is not exceeded. 

The program then checks to make sure that at this overpressure the 
percent saved does not also exceed the next higher desired percent saved. 

In other words if a table overpressure value yields a probability of coverage, 
hence mortality, which saves at least the desired percentage of population 
but does not save as many as the next higher desired percentage then that 
value of coverage probability is used to compute the number of survivors 
and the cost of shelters of the associated overpressure. 
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In case the overpressure found initially gives a value of probability 


of coverage outside the interval between the two mortality fractions, an 
interpolation is made on the radius of effect to reduce the overpressure so 
that the probability of coverage comes into this specified interval. It 
should be noted that in every case the overpressure finally found will be at 
least as great as that required to save the deBired percent exactly. The 
total numbers of survivors in the printout for a desired percentage will 
therefore always be greater than the desired percentage of the total population 
of that SLA. The cost of providing shelters of the design overpressure calcu¬ 
lated in the above manner is then found for the tract. * This procedure is 
repeated in the other tracts and the cost is summed over all tracts to find 
the total cost of providing the desired survival probability to the area being 
considered. 

This cost is then divided by the expected total number of survivors 
to obtain a cost-per-survivor value for each desired survival probability. 

A computer printout for a sample problem is reproduced in Table 
XXI. The printout provides an adequate description of the problem. The 

* 

In the present version of the program, it is assumed that shelters having 
a capacity equal to the tract population would be built. The program could 
be modified to accommodate some standard shelter capacities, or, two 
runs could be made. The first, with the real populations, could be used to 
estimate casualties; the second, with standard shelter capacities in place 
of tract population, cculd be used to estimate cost. 
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36 

cost table, COST (3), is based on a Guy B, Panero study; any other cost 
table of interest can be used in place of this one. 

The cost-per-survivor rows indicate that a minimum exists at about 
65 percent saved (60 percent desired) for COST (1), 74 percent for COST (2), 
and 87 percent for COST (3). This helps illustrate the importance of under¬ 
standing costs. However, the reader is cautioned against drawing general 
conclusions from this one printout; its primary purpose is to illustrate use 
of the program. 

3. 0 Uniform Shelter Design Overpressure 

This program permits calculation of some implications of a policy in 
which everyone is given space in shelters having Hie same design overpres¬ 
sure. The same alternative aiming points are allowed as in the uniform sur¬ 
vival probability program. Then, for the given attack variables, the computer 
calculates the probability of covering a geographic subarea (e. g. , SLA) with 
a particular overpressure. This probability is multiplied by the number of 
people present in the subarea, the process is repeated for all of the subareas, 
and the sum of the expected number of fatalities is calculated and printed for 
the chosen overpressure. The percent fatalities, the total cost of supplying 
everyone with sppcc ir. rheltrr of this design overpressure and the cost per 
survivor are also calculated and printed. The process is then repeated for 
each of the other overpressures supplied as input data. 

36 

Guy B. Panero, Inc., Shelter Configuration Factors , Prepared under 
Contract OCD-OS-62-108; April 15, 1963. 
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A sample problem is illustrated in Table XXII. The problem is ex¬ 
plained in the table heading. The point made in Par. 2. 0 about cost uncer¬ 
tainties can be made again with Table XXII. The overpressures associated 
with the minimum cost per survivor and total program costs are shown in 
Table XXIII. 


Table XXIII 

Effect of Cost Uncertainties in a Uniform Shelter Design 
Overpressure Sample Problem 


Cost function 

Overpressure (psi) 
for minimum cost 
per survivor 

Min. cost ($) 
per survivor 

Total cost 
($ millions) 

Expected % 
survivors 

COST (1) 

15 

204 

35 

57 

COST (2) 

20 

212 

42 

66 

COST (3) 

50 

251 

63 

84 


4. 0 Maximum Added Survivors per Dollar 

This program makes use of the same input data as in Pars. 2 and 3 
of this chapter. The probability of coverage is computed for each geo¬ 
graphic subarea (SLA) for each of a number of assumed shelter design over¬ 
pressures. Next, it is assumed that everyone in the geographic area (e. g., 
city) is given at least the protection of a shelter having some low value of 
design overpressure, referred to as the BASE overpressure. The expected 
number of survivors added by raising the shelter design overpressure to 
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each of the other overpressures chosen is then calculated for each SLA and 
divided by the cost of upgrading the protection to these chosen levels. The 
maximum value is found for each SLA and these maximum values are then 
ordered in decreasing order of expected added survivors per dollar. The 
value at the top of the list is printed along with identifying information. This 
value will be associated with a particular SLA and increase in shelter design 
overpressure. For example, suppose that upgrading the shelter in SLA 58 
from 1 to 3 psi yields the maximum expected added survivors per dollar. 

Then, SLA 58 must now be treated as if its BASE overpressure were 3 psi 
and the number of added survivors per dollar must be computed for raising 
SLA 53 shelters from 3 psi to each of the higher shelter design overpressures. 
The maximum of these values is found and replaces the old maximum value 
for SLA 58. The whole list is reordered, including this new value, its 
maximum value is printed, and the process is repeated until a complete 
printout is obtained of the order of spending that maximizes expected added 
survivors per dollar. The list may be used, for example, to decide on the 
shelter design overpressure for each SLA for a fixed budget. The printout 
contains an accumulated cost column. When the cost of supplying shelters 
of the original BASE overpressure is added to any particular value in this 
column, the total cost of the program is then obtained for the shelter system 
indicated to that point. Then, all that is needed is to enter this total cost 
list at the budget level of interest and proceed toward the beginning of the 
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printout, selecting the maximum shelter design overpressure found for each 
SLA. The list will also provide an order for construction that leads to the 
greatest number of survivors per dollar spent up to any particular time. 

A 3ample problem printout is shown in Table XXIV; the heading ex¬ 
plains the problem. The probability of coverage table is printed because it 
is useful for checking purposes. 
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CHAPTER V 


A MODEL FOR DEVELOPMENT OF PREFERRED MIXTURES 
OF EVACUATION AND SHELTER 

1.0 Introduction 

This chapter is concerned with the development of a mathema¬ 
tical model for use in attempts to delineate preferred mixtures of shelter 
and evacuation. Warning time is considered more explicitly and in more 
detail than in past work. 

Only one target area is considered in the model; the entire threat 
is considered to be that due to a nuclear attack on the area. It is assumed 
that the location of the aiming point is known and that one weapon is as¬ 
signed to the aiming point. The weapon parameters (yield, burst height, 
and delivery accuracy) and the distribution of possible times of detonation 
are input parameters to the model. 

It is assumed that shelters can be placed at any location. Shelter 
cost is assumed to increase with increasing hardness; hardness is as¬ 
sumed to be a measurable, reproducible quantity which indicates the vul¬ 
nerability of those in the shelter. The evacuation rate and vulnerability 
of ewcuees must be supplied as input parameters to the model. 
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2. 0 The Mathematical Model 


A path is defined over which it is planned to evacuate people from 
the target area. Along this path a distance coordinate s is defined. For 
each point along the path, there is a survival probability for an evacuating 
person who is at that point at the time of a detonation. This probability of 
surviving an attack while traveling is denoted P^,(s). The function P,p(s) 
depends upon delivery accuracy, weapon parameters, and the vulnerability 
of the evacuee. 

Along the evacuation path a distribution of shelters is assumed such 
that at any point an evacuee can stop and quickly enter a shelter. If a wea¬ 
pon detonation occurs after he has entered a shelter at point s, then his 
survival probability is denoted Pg(s) . This continuous distribution of 
shelters is hypothetical; the desired distribution of shelters and their hard 
ness are viewed as the unknowns in the problem. A continuous hypotheti¬ 
cal distribution of shelters is, however, a useful conceptual tool. 

Since a giv< n amount of protection becomes cheaper to attain as 
the distance from an aiming point increases, it is assumed that some peo¬ 
ple in the target area begin evacuation as soon as an alert is given. A 
typical evacuee undergoes motion described by a function of time s(t) 
along the evacuation route. Since he would be much safer in a shelter at 
the time of detonation than if he were traveling, he should stop and enter a 
nearby shelter at an appropriate time t . His survival probability for the 
attack is then 
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t 

P sur * S ° P Tt s(t) ) TU)dt * V^o’l I" ™ dt ' 

O t 


where T(t) is defined by the following three properties: 

pt+At 

(1) \ T(t ) dt is the probability of a detonation occurring between 

J t 

times t and t + At , 


T(t) dt = 1, and 


(3) * lm T(t) = 0 . 
t 


If P (s) , P (s), and T(t) are all considered to be known functions, 

J 1 

the maxima and minima of the survival probability with respect to the take- 
shelter time can be found by setting dP ^ I dt equal to zero and simplifying: 


P* <s ) s' (t » f T(t) dt - [P_(s ) - )] T(t ) -- 0 , 

So °^t l So io J o 

o 


where the primes denote differentiation, and s = s(t ) Thus, solving 
r o o 

Eq. (2) for t vields the time at which a person should stop traveling and 
o ' 

enter a shelter to make his survival probability a maximum. 

On the other hand, Eq. (2) can be viewed -•*£ a first-order linear 

ordinary differential equation in P (s ). Functions P satisfying Eq. (2) 

b O o 

for all s have the property that they make P constant, that is, 
o r sur 
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independent of the take-shelter time t . Such a shelter distribution has 

increasing safety in the direction of motion such that the added risk of being 

exposed long enough to move an additional short distance is always exactly 

offset by the safety of the shelter entered after the movement. The general 

solution of Eq. (2) for the function P can be written; 

o 


D 

" S 


(s ) 
o 



where Q = Pg(sJ - ^ 


o T(t Q ) P T <S Q )ds o 
5 i V(S 0 


(3a) 


(3b) 


and where t is the time of arrival at the point s , s. is the initial value 
o r o 1 

I I 

of s , and v(s q ) is the travel velocity at the point s q . 

The general solution of Eq. (2) given by Eqs. (3) can be divided into 


three cases depending on whether the value of P (s.) is less than, equal to, 

O 1 


or greater than a , 


w'here a 


r 

s. 

1 


T(t’) P (s )ds* 
_ o 1 o _o 

v(s ) 
o 


* 

Case 1, P^ (s ) < a . For some value of s , say s , the quantity Q 
Si o o 

* 

vanishes and tor s q > s^ it becomes negative. Thus, the only physically 

& 

allowable values of P 0 (s ) are for s < s (see Fig. 17). Following the 

So o — o 

discussion of Case 2 it will be shown that stopping in any shelter in a 
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distribution of shelters described in this case would result in lower survival 


probability than that associated with continuing the assumed motion given by 
s(t). 

f 00 

Case 2, P_(s.) = a . Let b be the least value of s such that \ T(t)dt - 0 
Si o J t 

o 

(if b is not finite, the limit as s -<x is to be substituted for s -»b in the 

o o 

following). As s^ approaches b, the quantity Q approaches zero. Since 

pOO 

the denominator of Eq (3a), \ T(t)dt, approaches zero at the same 

*o 

point, an indeterminate form is obtained for P (b). In evaluating the limit 
lim 6 

one finds that s —b P„(s ) = P_(b). Thus the solution given in this case 
o bo T 

can be physically realized for all values of s < b. (P„(s ) for s > b is 

o~ So o 

of no interest). Such a system of shelters has the property that stopping m 
any one of them gives the same survival probability as total evacuation by 
continuation of the assumed motion s(r), 

Note that P for Case 1 is always less than for Case 2 which pro- 

O 

vides the shelter distribution that has survival probability equal to that 
associated with evacuation. Thus, Case 1 shelters are less safe than eva¬ 
cuation 

Case 3 ; P (s )> o . In this case, the quantity Q is positive for ail values 

c x 

of 3 Since the denominator of Eq (3a) \ T(t) dt. approaches zero for 

° v t 

o 
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some value (finite or infinite) of t and hence for some value of s , say 

o o 

s^ = b, the solutions have a pole at S Q s b. Thus for s q greater than some 

v > r 

value, say s , P_(s ) is greater than unity and has no physical meaning. 

O D O 

The portions of the solution for which Pg < 1 represent shelter distri¬ 
butions that offer greater safety than an attempt at total evacuation; how¬ 
ever, notice that they exist physically for only limited distances from the 
target area. 

Figure 17 is a sketch showing the features of the functions described 
under Cases 1, 2, and 3 above. 


Case (3) solutions 



Fig. 17 

Sketch oi P_ and P_ for Various Constant Values of P 

T S sur 


3. 0 Preferred Evacuation-Shelter Policies for Given Shelter Systems 
The general solution of Eq. (2) for Pg(S Q ) has a usefulness aside 
from its direct application to he shelter system design problem. Suppose 
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a function P g <s q ) is given as a result of some other considerations (possibly 

it corresponds to an existing shelter system or was theoretically derived 

or. the basis of some other theory). It is then desired to find points s q for 

which Eq. (2) is satisfied and thus find the optimum locations for evacuees 

to take shelter. At a point where Eq. (2) is satisfied, the given function 

P (s ) must be tangent to one of the curves belonging to the general solu- 

tior. This is true because it must have a value and slope satisfying Eq. (2). 

Thus for a given T(t) and P r ^,(s), if a grid is constructed of the general 

solution curves of Eq. (2) and then some other P (s ) is plotted on the same 

ij o 

coordinates, the tangency points can be seen and thus (at least approximately) 
the points of solution of Eq. (2) can be read off the graph. Furthermore, in 
regions where the given function P Q (s ) crosses the general solution grid 
lines upwardly for increasing s q , a traveler will increase his safety by pro¬ 
ceeding farther. Similarly, in regions where P (s ) crosses the general 

o O 

solution grid lines downwardly for increasing s q , a traveler would decrease 
his safety by continued travel. 

4. 0 A Motivational Problem 

When people in the same area have instructions to do different things 
to increase their safety they will probably question whether they or some 
other group are increasing their safety most effectively. The self- 
preservation instinct provides a strong motivation to do that which is thought 
to be safest. Thus, people instructed to stay in a shelter may be tempted 
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to join with others who are traveling by just outside their door on their way 
to a shelter presumed to be safer because it is farther from the target area. 
On the other hand, those who are traveling may feel that they are much too 
vulnerable while traveling and therefore that they would be safer in one of 
the nearby shelters. Yet the civil defense plan may be designed to utilize 
the transportation facilities to the fullest leaving no transportation for those 
instructed to stay in ohelters. Conversely, the shelters along the evacua¬ 
tion path may have no room for those instructed to travel, shelters for the 
travelers having been provided farther from the target area where a given 
shelter hardness provides a greater survival probability. 

Though there would probably be some discontent in a situation such 
as the one des cribed no matter what shelter plan had been followed, the 
cause for the discontent can be eliminated, in principle at least, by using 
a shelter distribution with safety given by Cases (2) or (3) of the general 
solution of Eq. (2) for P (s ). If the shelters correspond to Case (2), the 

u O 

travelers can evacuate completely and those in the shelters will have the 
same survival probability as those that have evacuated. If the shelters 
correspond to Case (3), the travelers will eventually have to stop and enter 
a shelter, but those who stop and enter a shelter first will have the same 
survival probability as those who go farther. 

One difficulty soon becomes evident in realizing shelters that corres¬ 
pond to Cases (2) or (3); these theoretical shelter distributions are, in 
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general, dependent on the initial position, s . of the traveler To circum¬ 
vent this problem, a mean value of s. for a traveling group can be used, or 
a special form of T(t) can be used which removes the dependence on s . 

Such a form for T(t) is presented in the section in which the function T(t) 
is discussed. In spite of this difficulty the theoretical distributions given by 
Cases (2) and (3) are of considerable conceptual significance 

5. 0 Evacuation-Shelter Policy for a Fixed Budget 

Shelter-building funds might be appropriated and distributed on the 
basis of population density. If so, a given population density would imply a 
particular shelter hardness, represented by shelter design overpressure. 

To locate shelters for small groups having a fixed budget, it is only neces¬ 
sary to construct the function P (s ) as a function of distance along the 

o O 

evacuation path. This function is inserted in Eq. (2) which is then solved 
for the maximum of total survival probability. Mechanically, this can be 
done by plotting P (s ) on the same graph with members of the set of func- 

o O 

tions that represent the genera] solution of Eq (2) and then looking for tan¬ 
gent points. If there is more than one tangent point representing a maximum, 
a shelter should be built where it offers the greatest survival probability. 
Figure 18 illustrates the technique 

6. 0 Some Comments On and Examples Of The Function T(t) 

An important property of the curves belonging to the general solu¬ 
tion of Eq. (2) , given by Eqs. (3) , should be noticed. In regions where 
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target area Distance from aiming point 


Fig. 18 

Illustration of Technique for Finding the 
Optimum Location for a Shelter of a Given Cost 

poo 

T(t ) = 0, but \ T(t) dt > 0 , the function P c (s ) = a constant, 
o J bo 

o 

Thus, the initial portion of the curves given by Eq. (3) is constant 

ar.d equal to P r ,(s.) over the distance traveled by the person in question 
o i 

during any hiitial time in which T(t) = 0 . 

In the examples which follow it will be assumed that the motion of 
the traveler is rji.en by s{t) = +• vt, where v is a constant. Specific 
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forms of the function T(t) will be chosen; they a*e not necessarily thought 
to be realistic, only instructive. 

Example (1) 


Let T{t) =0, t < ♦ 


= ( 1 /t) exp ^-(t-tj)/^ , t > tj 

For t < t, , that is for s. < s < s. + vt. , P_ (s ) = P e (s.) 

1 i-oi ISo Si 


For t > tj , 


P e (s ) = exp (s / vt) | P_ (s.) exp (-s. / vr - t. / r) 
So o {. S i r i 1 


- U / vt) 


\ Pm(s ) exp (-s / vt) ds \ 

J „ To ^ o o I 


s. 

1 


This can be rewritten 


P (s ) - exp (s / vt) { C -(1 /vt) \ P (s ) exo(-s /vT)ds ] , 

So F o l J 0 To o o/ 


(4) 


where C = P_ (s.) exp (-s./vr - t /t) + (1 / vr) 

O 1 1 i 


S . 

I 1 
J 0 


P_(s )exp(-s /vT)ds 
Too o 


The family of curves given by Eq. (4) is the general solution of Eq. (2) and 

is independent of s . The constant of integration is related to ar shown, 

but the same values of C are covered independent of s since P_(s.) is an 

arbitrary constant. It should be remembered, however, that the solutions 

given by Eq. (4) hold only for t > t ; this implies that s > +■ vt . 

I o ** i i 
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This is simply a reminder that the first portion (s^ < s > vt^) of each 
curve is a constant and then the curve should join one of those given by 
Eq. (4). Thus, if the family of curves given by Eq. (4) is plotted, an actual 
solution of Eq. (2) for a given s. is formed by joining a horizontal straight 
line (P Q = P c (s.)) drawn from s. to s. + vt , to the proper curve given 
by Eq. (4). Figure 19 illustrates this construction. 



target area Distance from aiming point 

Fig. 19 


Construction of Solutions of Eq. (2) from 
the Family of Curves Given by Eq. (4) 

The particular form of T(t) discussed in this example has the inter 
esting property that the general solution of Eq. (2) is independent of s. 






except for the (probably) small portion between s and *■ vt If t -0. 
the general solution is completely independent of s . 

Example (2) 

Let TU) =0, t < t t > t 

= 1 /(t 9 - t 1 >, t 1 < t < t 2 

As in Example (1), P (s ) = P (s ) for t < t (or equivalently for 

b o b 1 1 

s. < s < s. + vt,) . The general solution of Eq. (2) for t, < t < t_ can be 
l “ o l 1 b 1 ~ 2 

written (from Eq (3) ). 


P Q (s ) 

S o 


2 1 


t - - (s -s.) ) S 
2 V o l ! 


!’(«•) - 77 
b l (t 


~— r 

-tjv J 


2 1 


s.-t-vt, 

l 1 


P (s ) fis ) , (5) 

I o of 


where s. f vt, < s < s. + vt,, The function P, (s ) is physically meaning- 

l 1 — o i 2 bo 

less for t > t„ (or s > s vt, ) since there is no longer anv threat of 
2 o l 2 ^ 

attack hence P ( , (s ) will not be discussed further for that domain 
b o 

Ail solutions are constant for s < s < s. \t. Then the solution 

l — o l 1 

functions P (s ) begin to increase, taking on the values given by Eq (5) 

Case (2) solutions intersect P.„ (s ) at s = s + vt„ and Case (3) solu- 

1 o o l 2 

tions are singular at that point becoming greater than one for some s^ < 

* vt (see Eig 2Ui I'nlikc the functions P e (s ) of Example (1), the solu- 
2 bo 

tions P.(s ) in this example are dependent on s 
b o l 

A computer program to solve Eq (2) has been developed and applied 
to some practical problems, a sample of the results obtained is provided in 


Table XXV and in I :g 21 
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Survival probability 


(5) 


/ 

/ 



target area Distance from aiming point 

Fig. 20 

Illustration of the General Solution of Eq. (2) for the 
Function T(t) as defined in Example (2) 






Typical Computer Output for Program * b ; M * r 

Mixtures for Constant Overall Survive* Pr^t.ann ;*v 

PROGRAM GOSH 

K. D. GRANZOW, OIKEWCOO CORP 
MAY 15* 1964 

OS « 0,500 MILES NOS * 20 

CEP « 2.00 MILES YIELD * 10.00 MEGATON 

VELOCITY * 3.00 MILES/HR 

SI • 0. MILES 

ATTACK PROBABILITY ZERO FOR T LESS THAN 0.2 HOURS 

ATTACK PROBABILITY EXPONENTIAL WITH TIME CONSTANT 1.0 HOURS 


HARONESS LIST IP 

SI) - 

1, 2. 5, 10. 20, 50, 100, 200 

-CONSTANT SURVIVAL PROBABILITY- 

S 

PT(S) 

PS<S1(CASE 2)/ PS< S) ✓ 

OVERPRESSURE OVERPRESSURE 

PSISI) * 

HARONESS 


0.19- 

0.32 0.46 0.59 0.73 0.86 

0. 

0.00 

0.19- 

20.0 

0.32 0.46 0.59 0.73 0.86 

30.0 40.0 75.0 150.0 500.0 

0.5 

0.00 

0.19 

20.0 

0.32 0.46 0.59 0.73 0.86 

30.0 40.0 75.0 100.0 500.0 

1.0 

0.00 

0.22 

20.0 

0.37 0.52 0.68 0.83 0.99 

30.0 50.0 75.0 200.05000.0 

1.5 

0.00 

0.25 

20.0 

0.44 0.62 0.80 0.98 1.17 

30.0 50.0 100.05000.06000.0 

2.0 

0.01 

0.30 

15.0 

0.51 0.73 0.95 1.16 1.38 

30.0 75.0 500.06000.06000.0 

2.5 

0.02 

0.35 

15.0 

0.61 0.86 l.ll 1.37 1.62 

25.0 100.06000.06000.06000.0 

3.0 

0.04 

0.41 

15.0 

0.71 1.01 1.31 1.61 1.91 

30.06000.06000.06000.06000.0 

3.5 

0.07 

0.47 

15.0 

0.83 1.18 1.54 1.90 2.25 

30.06000.06000.06000.06000.0 

4.0 

0.11 

0.54 

15.0 

0.96 1.38 1.80 2.22 2.64 

100.06000.06000.06000.06000.0 

4.5 

0.18 

0.62 

10.0 

1.11 1.61 2.10 2.60 3.10 

6000.06000.06000.06000.06000.0 

5.0 

0.26 

0.69 

10.0 

1.27 1.86 2.45 3.03 3.62 

6000.06000.06000.06000.06000.0 

5.5 

0.36 

0.76 

8.0 

1.45 2.14 2.83 3.53 4.22 

6000.06000.06000.06000.06000.0 

6.0 

0.48 

0.62 

8.0 

1.64 2.45 3.27 4.09 4.91 

6000.06000.06000.06000.06000.0 

6.5 

0.59 

0.87 

8.0 

1.84 2.80 3.77 4.73 5.70 

6000.06000.06000.06000.06000.0 

7.0 

0.70 

0.91 

8.0 

2.05 3.19 4.34 5.40 6.62 

6000.06000.06000.06000.06000.0 

7.5 

0.79 

0.94 

8.0 

2.29 3.64 4.99 fc.33 7.68 

6000.06000.06000.06000.06000.0 

8.0 

0.86 

0.96 

8.0 

2.56 4.15 5.74 7.33 8.93 

6000.06000.06000.06000.06000.0 

8.5 

0.92 

0.98 

8.0 

2.86 4.74 6.62 8.50 10.38 

6000.06000.06000.06000.G6000.0 

9.0 

0.95 

0.98 

8.0 

3.21 5.43 7.65 9.87 12.10 

6000.06000.06000.06000.06000.0 

9.5 

0.98 

0.99 

6.0 

3.61 6.24 8.86 11.49 14.12 

6000.06000.06000.06000.06000.0 

10.0 

0.99 

0.99 

5.0 

4.09 7.19 10.29 13.40 16.50 

6000.06000.06000.06000.06000.0 


500 
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Survival _>ro*:a-.ilit.y Overpressure (psi) 


1000 


86 


Required Shelter 



Distance from aiming point (miles) 

Note: Numbers on curves are overall 
survival probabilities. 



Fig. 21 

Constant Survival Probability Curves tor 
Mixtures of Movement and Shelter 
-128- 


0, t<0. 2 hour 
negative 
exponential with 
time constant of 
1 hour for t>0. 2 
hour 
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ABSTRACT 


This report describes an effort to find preferred mixtures of 
movement and shelter as Civil Defense responses to the threat of nuclear 
war. Two approaches were followed 

1. Mixtures of movement ai d shelter were studied in three 
steps. These consisted of 

a. Postulation of alternative movement and shelter policies, 

b. Development of movement and shelter plans based on 
these policies, 

e. Evaluation of plans developed in (b) against the range of 
attack conditions considered reasonable. 

2. A mathematical model was constructed to provide a 
vehicle for sensitivity analyses. 

A technique for planning large -scale strategic movements was 
developed and applied to several particular places. The technique is 
believed to be developed sufficiently to provide a basis for planning a 
first - generation strategic movement capability for the U. S. 

Two computer programs were developed as tools for evaluating 
strategic movement against particular attacks and for evaluating var¬ 
ious trans-attack responses to large-scale movements interrupted by 
war. 

Blast shelter planning programs are also reviewed and developed 
furthe Evaluation techniques are already available. 

The mathematical model aoproach ended with the development 
of a computer program for finding the shelter location and hardness re¬ 
quired to maximize overall survival probability for various warning 
i.ime probabil iy density functions and attack and movement assumptions. 
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